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In this paper, as a personal review, we suppose a possible extension of Gibbs ensemble theory so 
that it can provide a reasonable description to phase transitions and spontaneous symmetry break- 
ing. The extension is founded on three hypotheses, and can be regarded as a microscopic edition of 
the Landau phenomenological theory of phase transitions. Within its framework, the stable state 
of a system is determined by the evolution of order parameter with temperature according to such 
a principle that the entropy of the system will reach its minimum in this state. The evolution of 
order parameter can cause change in representation of the system Hamiltonian; different phases 
will realize different representations, respectively; a phase transition amounts to a representation 
transformation. Physically, it turns out that phase transitions originate from the automatic interfer- 
ence among matter waves as temperature is cooled down. Typical quantum many-body systems are 
studied with this extended ensemble theory. We regain the Bardeen-Cooper-Schrieffer solution for 
the weak-coupling superconductivity, and prove that it is stable. We find that negative-temperature 
and laser phases arise from the same mechanism as phase transitions, and that they are instable. 
For the ideal Bose gas, we demonstrate that it will produce Bose-Einstein condensation (BEG) in 
the thermodynamic limit, which confirms exactly Einstein's deep physical insight. In contrast, there 
is no BEG either within the phonon gas in a black body or within the ideal photon gas in a solid 
body. We prove that it is not admissible to quantize Dirac field by using Bose-Einstein statistics. 
We show that a structural phase transition belongs physically to the BEG happening in configu- 
ration space, and that a double-well anharmonic system will undergo a structural phase transition 
at a finite temperature. For the 0(A'^)-symmetric vector model, we demonstrate that it will yield 
spontaneous symmetry breaking and produce Goldstone bosons; and if it is coupled with a gauge 
field, the gauge field will obtain a mass (Higgs mechanism). Also, we show that an interacting Bose 
gas is stable only if the interaction is repulsive. For the weak interaction case, we find that the BEG 
is a "A" -transition and its transition temperature can be lowered by the repulsive interaction. In 
connection with liquid ''He, it is found that the specific heat at constant pressure Cp will show a T"^ 
law at low temperatures, which is in agreement with the experiment. If the system is further cooled 
down, the theory predicts that Cp will vanish linearly as T — > 0, which is anticipating experimental 
verifications. 

PACS numbers: 05.30.Ch, 05.30.-d, 64.60.-i 
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1. INTRODUCTION 
1.1. The Problem of Gibbs Ensemble Theory 

Since Gibbs [J| established the ensemble theory of sta- 
tistical mechanics, there has arisen the problem as to 
whether it has the ability to describe phase transitions 
0. 

Let us investigate the problem through the famous ex- 
ample, i.e., the conventional low-Tc superconductivity, 
which can be described by the so-called Bardeen-Cooper- 
Schrieffer (BCS) Hamiltonian 

H{c) ^ £(k)(c[^CkT + c^^kiC-ki) 



'Zl^k'T^-k'i'^-kiCkT, 

k,k' 



(1.1) 



where e(k) denotes the energy relative to Fermi level, 
g > the coupling strength, Ck| and Ckj (cj^| and cjj.|) 
the destruction (creation) operators for the electrons with 
up and down spins, respectively. Here, the notation H{c) 
expresses that the Hamiltonian is a function of Ck| , Ckx , 



Ck| , and c|j.| . 



For BCS superconductivity, what is the most impor- 
tant physically is that the Hamiltonian H{c) will remain 
invariant under the transformation. 



G{d,c)H{c)G\'d,c)^H{c), 



(1.2) 



where 



G{d, c) = e-'''^>'("M^>'T+4i=ki)_ ^ g [Q^ 27r). (1.3) 

This invariance is known as gauge symmetry. As a direct 
consequence of this symmetry, one has 

(c^kiCkt) = Tr(c_kiCkTP(iJ(c))) = 0, (1.4) 

where (■ • • ) denotes the Gibbs ensemble average with re- 
spect to the Hamiltonian -ff(c), and 



P{H) 



(1.5) 



is the statistical density operator where (3 = 1/ (ksT) 
with ks and T being the Boltzmann constant and tem- 
perature, respectively. To prove Eq. H1.4|) . it is sufficient 



to heed that 



C-kiCkiG^ (^,c) = -c_kiCkT- (1-6) 



Eq. (|1.4I) indicates that the electron-pair amplitude 
(c_kiCk|) can never become nonzero according to Gibbs 
ensemble theory, in other words, there can not happen 
superconductivity for the BCS Hamiltonian of Eq. Hl.l() 
within Gibbs ensemble theory no matter how low the 
temperature is. This proves definitely that Gibbs ensem- 
ble theory has no ability to describe phase transitions. 

Even so, it is still hoped to do, as least as possible, some 
modification to Gibbs theory so that it could describe 
phase transitions. Up to now, there are two kinds of 
modification, they both arise from the theoretical studies 
on the so-called Ising model 6] , and become the popular 
and predominant believes about phase transitions in the 
current world of statistical physics. 

The first is to introduce the thermod yna mic limit in 
the end of calculation 0, i, H US [H E llllll , that is. 



(c-kiCkt; 



lim Tr(c_kiCkTP(-H^(c))), 

V— ► + 00 



(1.7) 



where the limy_>4.oo means to take the thermodynamic 
limit {V +00, N +oo, and N/V = constant) after 
performing the trace. Because 



Tr(c_kiCkTP(i/(c))) = 0, 



(1.8) 



irrespective of the magnitudes of V, N, and N/V {H{c) 
is invariant under gauge transformation, and that is true 
independently of V^, N, and N/V, namely, the gauge 
symmetry has nothing to do with V , N, and N/V.), one 
has 



lim Tr(c_ 



ki 



(1.9) 



Therefore, this kind of modification can not make Gibbs 
ensemble theory have the ability to describe phase tran- 
sitions. Indeed, the thermodynamic limit is important 
for calculating the statistical average of observable, but 
itself alone is insufficient to act as the physical criterion 
for phase transitions. 

The second is to introduce both the thermodynamic 
limit and an auxiliary external field 



the tnermodyn 



— lim lim Tr 

0^0 V^+oo 



-kiCkTp(i?(c) + ^((/)C_kiCkT + ^^cj^T^-ki)) 



(1.10) 



where, as pointed out by Huang [17| . the limit — > must be taken after the thermodynamic limit V +oo. 



3 



The purpose of introducing an auxihary external field (j) 
is to break the gauge symmetry of the system, as can be 
easily seen from the term, 

^(0c_kiCkT + </'^4^c^_ki)- (1-11) 

k 

Since the gauge symmetry has been broken with regard 
to the entire Hamiltonian, 

H{c) + ^(0c_kiCkT + -^^^Tclkj,)' (1-12) 

k 

in so far as the auxiliary external field (j) is infinitesi- 
mal but nonzero, the trace and the limit V — > +oo of 
Ea. (|l.lU|l become nonzero. Further, if the limit — s- 
does not tend to zero, the electron-pair amplitude 
(c_kxCk|) will get nonzero. One thus concludes that there 
appear Cooper pairs in the system, and that the sys- 
tem has gone into the superconducting phase. However, 
this conclusion can not hold in physics, for the modifi- 
cation suffers the serious problem: The limit procedure 
employed by Eq. H1.10() is physically equivalent to the sce- 
nario that an infinitesimal auxiliary external field is first 
added to the system so as to induce the symmetry of the 
system to break down, and taken off from the system 
finally. But, in the actual situation, the inducement of 
such an external field is unnecessary, the symmetry it- 
self breaks down spontaneously and does not need any 
help of external force. Anyhow, this kind of modifica- 
tion contains unphysical operations, and thus can not be 
accepted in principle. 

Here, it is also significant and worth while to 
give a brief discussion to the two-dimensional Ising 
model because it can be solved exactly within the two 
modifications, 

Hl = -jY,SrSj, (1.13) 

where = ±1 denotes the spin variable on site i, the 
symbol {ij) means that the sites i and j are the nearest 
neighbors, and J > is the exchange coupling between 
a nearest-neighbor pair of spins. 

As is well known, Onsager proved rigorously that 
the specific heat of the Ising Hamiltonian H\ is singular 
at T = Tc > within the first modification of Gibbs 
ensemble theory. It hints that the system might undergo 
a paramagnetic-ferromagnetic phase transition at in 
the thermodynamic limit. However, just as Eqs. (|1. Ti - 
ll. 9|l . the magnetization to, 

m = lim Ty{sw{Hi)) = 0, (1.14) 

is always equal to zero due to the fact that the parity 
symmetry of Hi with respect to = 1 and Si = —1 has 
nothing to do with the magnitude of N where N repre- 
sents the total number of sites. Eq. H1.14|l proves defi- 
nitely that there is no paramagnetic-ferromagnetic phase 



transition at any temperature within the first modifica- 
tion of Gibbs ensemble theory even if the specific heat is 
singular. 

The above proof also shows that it can not be solved 
within the framework of the first modification whether 
or not the onset of the singularity at T = Tc manifests a 
phase transition. In order to justify that the phenomenon 
occurring at T = Tc is a phase transition, Yang 
went beyond the first modification and turned to the sec- 
ond one, which was, in fact, suggested originally by Ising 
himself 6] in studying the one-dimensional Ising model. 
Wonderfully, Yang succeeded in proving rigorously that 

TO = lim lim T!t\ SipiHi — By s.;) | 

i 

= 1^^™' (1.15) 
1^ nonzero, 1 < Ic, 

where B stands for an external magnetic field. Math- 
ematically, it seems reasonable to say from Eq. ()1.15|) 
that a paramagnetic-ferromagnetic phase transition will 
occur at T = Tc. Nevertheless, as pointed out above, 
Eq. (|1.15|) imports unphysical operations from the in- 
troduction of the external field i?, and this leads to the 
result that the phase transition will not occur automati- 
cally but has to be driven and decided from outside the 
system, which directly contradicts the basic experimental 
fact that any phase transition occurs itself spontaneously. 
Therefore, the second modification, or rather Ising's cri- 
terion, can not be used as a physical criterion to justify 
whether there exists a paramagnetic-ferromagnetic phase 
transition at T — Tc. 

It should be stressed again that it is the two rigorous 
works of Refs. 

and that establish the two kinds 
of modification to the Gibbs ensemble theory. 

On all accounts, the two kinds of modification must 
both be discarded, and it is necessary to modify or ex- 
tend Gibbs ensemble theory anew. That is just the main 
purpose of this paper. 

Recently, Gibbs ensemble theory has been extended 
by C. Tsallis and his followers with the conception of 
nonextensive entropy 0, IT^ I2(tI | , which is being under 
controversy. In this paper, as a personal review, we would 
like to suppose another possible extension of the Gibbs 
ensemble theory so that it can provide a reasonable de- 
scription to phase transitions and spontaneous symmetry 
breaking (SSB). 



1.2. Landau Phenomenological Theory 

Also, there is a macroscopic theory for phase transi- 
tions, i.e., Landau phenomenological theory [2ll |. 

Phenomenological as it is. Landau theory succeeds in 
providing a unified picture for all the second-kind phase 
transitions, e.g., superconductivity, superfluidity, mag- 
netism, and structural phase transitions. Physically, 
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Landau picture consists of two basic notions: order pa- 
rameter and variational principle. Order parameter de- 
scribes the degree of order of the system, it is zero in 
the disordered phase, and nonzero in the ordered one. 
Variational principle yields the equation of motion of or- 
der parameter and controls the evolution of order pa- 
rameter with temperature, a system must arrive at the 
minimum of Helmholtz or Gibbs free energy in its sta- 
ble state. According to this picture, a system will evolve 
with the evolution of its order parameters, and produce 
phase transitions spontaneously, i.e., without any help or 
drive from outside the system. In short. Order param- 
eter and variational principle are two characteristics of 
Landau theory. 

On the other hand, Landau pointed out that a phase 
transition of the second kind reflects physically the 
change in symmetry of the system; the disordered phase 
has a higher symmetry than the ordered one. He found 
that this change in symmetry can be described through 
the representations of the symmetry group, a phase tran- 
sition corresponds to a representation transformation, 
different phases will realize different representations, re- 
spectively. Accordingly, Landau established the relation 
between the order parameter and the representation. Us- 
ing this relation, it can be easily determined from the 
variational principle which representation will be realized 
at a certain temperature. Representation transformation 
and spontaneous symmetry breaking are the other two 
characteristics of Landau theory. 

To sum up, a system will select and realize automat- 
ically different representations at different temperatures 
according to the variational principle of order parameter. 
That is the mechanism for phase transitions discovered 
by Landau theory, it crystallizes Landau's ideas about 
phase transitions: order parameter, variational principle, 
representation transformation, and spontaneous symme- 
try breaking. Theoretically, this mechanism gives a rea- 
sonable answer to the problem why a phase transition or 
SSB can happen at a certain temperature. In applica- 
tions, it also agrees quite well with various phase transi- 
tions, its effectiveness and universality are well known to 
condensed matter physicists. In a word. Landau's ideas 
on phase transitions are successful and of great value in 
physics! 

Regrettably, it is impossible to deduce Landau theory 
from Gibbs ensemble theory because there is no varia- 
tional principle of order parameter within the latter. This 
demonstrates again the deficiency of Gibbs ensemble the- 
ory. In fact, this deficiency was pointed out by Born and 
Fucks 12^] early in 1938. In 1937, Mayer ^] published 
his famous paper on gas-liquid transition, which is based 
directly on Gibbs ensemble theory. This work was seri- 
ously doubted and questioned by Born and Fucks [2^ : 
"How can the gas molecules 'know' when they have to 
coagulate to form a liquid and solid?" Obviously, owing 
to lack of Landau mechanism, it is hard for Gibbs en- 
semble theory to answer the question posed by Born and 
Fucks. 



Gibbs ensemble theory should be extended to incorpo- 
rate the Landau's ideas so as to describe phase transitions 
and answer the question posed by Born and Fucks. That 
is another intention of this paper. 

2. THE EXTENDED ENSEMBLE THEORY 

In order to describe phase transitions and spontaneous 
symmetry breaking, we shall, as a personal review, ex- 
tend Gibbs ensemble theory with three hypotheses, which 
can be stated, taking the BCS superconductivity as an 
instance, as follows. 

(1) . The system Hamiltonian is represented by 
ff'(0,c), 

H'{(j), c) = e'-°('^'")i?(c)e-'-°('^'"\ (2.1) 

where 

D{<P, c) = ^(0kC-kiCkT + <4^clki)- (2.2) 

k 

Here is an internal field of the system, it is also the 
order parameter for BCS superconductivity. Evidently, 
H'{(j),c) is not invariant under the gauge transformation 
of G{'d, c) as long as 7^ 0, namely, the gauge symmetry 
will be broken for H'{4>,c) if </) ^ 0. From now on, we 
shall call D{4>,c) the phase-transition operator, for sake 
of convenience. 

(2) . The statistical average of an observable F{c) is 
defined as 

.F(0, (3) = (F(c)) = Tt{F{c)p{H'{4,, c))) . (2.3) 

We remark that the average is now a function of both 
temperature and order parameter, which is expressed ex- 
plicitly by the two arguments of J^{(f>, (3). 

(3) . The order parameter (j) is determined by the min- 
imum of the entropy of the system, 

5S = 0, (2.4) 
A5 > 0, (2.5) 

where S denotes the entropy, 

5(0, /3) = {-Hp{H{c)))) 

= -TT{\n{p{H{c)))p{H'{^,c))). (2.6) 

As a function of the order parameter, the entropy controls 
the evolution of state of the system with temperature. 

The first hypothesis means that the system Hamilto- 
nian can take different representations at different tem- 
peratures, e.g., it can take the symmetric representation 
{(j) = 0) at a high temperature, and the asymmetric repre- 
sentation {(j) ^ 0) at a low temperature. Which represen- 
tation it will take is determined by the third hypothesis. 
After the representation is so determined, the statistical 
average can be calculated with respect to this represen- 
tation, as stated in the second hypothesis. 
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Obviously, the extended theory will reduce to the orig- 
inal one if = 0, that is to say, the original theory holds 
only for the normal phase of the system. The broken- 
symmetry phase will be described by the extended the- 
ory. 

Within the framework of the extended ensemble the- 
ory, a system will realize different representations of the 
same system Hamiltonian at different temperatures ac- 
cording to the principle of least entropy with respect to 
order parameter, this mechanism for phase transitions 
is, in spirit, the same as that given by Landau phe- 
nomenological theory; hence, it can be said that Lan- 
dau mechanism has been incorporated into the extended 
ensemble theory. Also, because the order parameter is 
an internal field of the system itself, symmetry break- 
ing will occur in a completely spontaneous way rather 
than forced by an external field, there is no unphysical 
operation within the extended theory. In a word, the ex- 
tended ensemble theory is conceptionally in accordance 
with Landau's ideas on phase transitions: order parame- 
ter, variational principle, representation transformation, 
and spontaneous symmetry breaking. 

In Sec. ini we shall show further that phase transitions 
originate physically from the wave nature of matter. 

Now, let us establish the relationship between the ex- 
tended ensemble theory and thermodynamics, it can be 
implemented as follows. 

First, the internal energy U can be obtained from Eq. 
(|2.3|) . replacing F(c) by the Hamiltonian H{c). It is a 
function of temperature T and volume V, that is, 

U = UiT,V). (2.7) 

From the internal energy U{T,V), the specific heat at 
constant volume, Cy, can be calculated through 

Cv^CviT,V)=(^^^^. (2.8) 

Then, we can obtain the thermodynamical entropy Sth, 

dT' 

Sth = SthiT, V) = Cv{T', V) — . (2.9) 

It is worth paying attention to the significant difference 
between the thermodynamical entropy Sth and the sta- 
tistical entropy S defining in Eq. (|2.6|) : S measures the 
degree of order of the system, whereas Sth measures the 
amount of heat absorbed or rejected by the system, i.e., 
dQ = TdSth] besides, S contains the basic information 
of the system, it controls the evolution of state of the 
system, as a consequence, Sth derives itself from S. 
At last, we arrive at 

F = F{T, V) - U{T, V) - TSthiT, V), (2.10) 

where F is the Helmholtz free energy of the system, with 
T and V as its natural variables. As a thermodynamic 
potential of the system, F can generate other thermody- 
namic quantities. By such a way, the extended ensemble 
theory connects with thermodynamics. 



In the end of this section, let us explain why we choose 
the BCS model rather than the Ising model as the start- 
ing instance to develop the theory for phase transitions 
and spontaneous symmetry breaking. 

As is well known, Ising model Hi is a classical discrete 
Hamiltonian. Simple as it is mathematically, it is, how- 
ever, ill defined from the point of view of physics. On 
one hand, a classical Hamiltonian, H = H{q,p) where q 
and p denote the generalized coordinates and momenta 
of the system, must be continuous in phase space, that is 
to say, H{q,p) is a continuous function of q and p; there 
does not exist any classical system which corresponds to a 
discrete Hamiltonian; therefore. Hi does not correspond 
to a classical system. On the other hand, the discrete 
quantity ±1 is not a quantum operator but a pure 
integer variable, therefore. Hi does not correspond to a 
quantum system, either. In a word, Ising model can not 
correspond to any physical system; its solution, whether 
rigorous or not, has no real physical relevance [21, which 
was, in fact, realized and pointed out as the "Ising dis- 
ease" by people early in 1940's and 1950's 

To understand that point further, let us first suppose 
that Ising model could represent a real physical system, 
and then analyze the statistical properties of this sys- 
tem in detail. First, if the system is not placed in an 
external magnetic field, then, according to the Onsager's 
solution 0, [13 , its specific heat will diverge logarithmi- 
cally at T = Tc > 0; however, as indicated by Eq. H1.14|l . 
there appears no spontaneous magnetization bellow Tc, 
namely, the state of the system will remain nonmagnetic 
at any temperature, whether T > Tc or T < Tc- There 
is no paramagnetic-ferromagnetic phase transition in the 
Ising system if it is not placed in an external magnetic 
field, though its specific heat is singular. That is the first 
characteristic behavior of the Ising system. Secondly, if 
the system is placed in an external magnetic field, then, 
according to the Yang's solution there will appear 
an induced magnetization bellow Tc- That is to say, the 
Ising system will transform from a paramagnetic state at 
T > Tc into a ferromagnetic state at T < Tc when it is 
placed in an external magnetic field, here, the transfor- 
mation is an induced transformation other than a sponta- 
neous phase transition. That is the second characteristic 
behavior of the Ising system. In sum, an Ising system 
will exhibit two distinct behaviors, respectively, depend- 
ing on whether the system is placed in an external mag- 
netic field or not. It will exhibit the first characteristic 
behavior in the case without the application of any ex- 
ternal magnetic field and the second one in the case with 
the application of an external magnetic field. The two 
behaviors should be checked experimentally with respect 
to the same physical system. However, such a physical 
system that can simultaneously show both the two char- 
acteristic behaviors of the Ising model has never been 
observed and reported experimentally. In other words, 
there does not exist the supposed Ising system in nature 
at all! That is not surprising, it just reflects the ill defi- 
nition of the Ising model. Here, it should be pointed out 
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that the solutions of Onsager and Yang are not identical 
but essentially different from each other, whether in the 
sense of mathematics or in the sense of physics: In the 
limit B ^ 0, the Yang's solution can not reduce to the 
Onsager 's solution, as can be easily seen by comparing 
Eq. (|1.14(l and Eq. H1.15(l : they represent the statisti- 
cal properties of the Ising system with and without the 
interaction of an external magnetic field, respectively. 

Needless to say, an Ising system is completely distin- 
guishable from a real physical system of phase transi- 
tion. For the latter, the transition happens itself spon- 
taneously, needing no inducement of a corresponding ex- 
ternal field. Evidently, this nature of phase transitions 
is neither the same as the first characteristic behavior 
of the Ising system, i.e., a singular specific heat but no 
phase transition, nor the same as the second character- 
istic behavior of the Ising system, i.e., a transformation 
induced but not spontaneous. There are clear and defi- 
nite discrepancies between the statistical behavior of a 
real physical system of phase transition and those of 
the Ising system. Any phase transition can not be de- 
scribed by Ising model. Of course, if one disregards 
physically the ill definition of Ising model, the distinc- 
tion between the solutions of Onsager and Yang, and 
the discrepancies between the statistical behavior of a 
real physical system of phase transition and those of 
the Ising system, he can simulate some phase transi- 
tions and critical phenomena quite well with Ising model 

[11 111 m H El El US El il . 

Theoretically, Ising model is ill defined; experimen- 
tally, Ising system does not exist in nature at all. That 
is the reason why we do not choose it as the starting 
instance to develop the theory of phase transitions and 
SSB. As for the BCS model (aHlll, quite the contrary, 
it is of high value in theory because it contributes a new 
quantum concept, i.e.. Cooper pairs, which is profound 



where 

/(^) = (3.5) 

is the Fermi distribution function. 

From Eqs. ()3.4|l and (|2.4() . it follows that 

sin(26'k) = 0, (3.6) 



and rather hard to understand theoretically but yet veri- 
fied exactly by the quantized flux experiments. Actually, 
it is the BCS theory that began our microscopic under- 
standing of phase transitions. 

Although Ising's criterion for phase transitions is prob- 
lematic in physics, it is much suggestive in mathematics, 
indeed. To some degree, it can be said that the present 
modification is just equivalent to substituting the exter- 
nal field proposed by Ising with an internal field deter- 
mined by the system itself, as can be seen from expanding 
the right-hand side of Eq. H2.1|) to the linear term of 4> 
and then comparing it with Eq. H1.12|l . 



3. APPLICATION TO THE IDEAL FERMI GAS 

Let us first apply the extended ensemble theory to the 
ideal Fermi gas. It is the simplest case, and we shall see 
that this case is exactly solvable. 

For the ideal Fermi gas, there is no interaction (g — 0), 
Eq. reduces to 

H{c) = ^e(k)(c[^CkT +c^_kiC-ki), (3.1) 

k 

where e(k) = ?i^k^/ (2m) — /i with m and ^ being the 
mass of the fermions and the chemical potential of the 
system, respectively. 

To facilitate the calculation of the entropy, we refor- 
mulated Eq. H2.6|l as 

5(^,/3) = -Tr(ln(p(i?(e-'^^(^^=)ce'^^(^^=')))p(i/(c))) . 

(3.2) 

By use of Eq. (|2.2|l . we find 



(3.3a) 
(3.3b) 

(3.4) 



which is the equation of order parameter. Obviously, this 
equation has the solutions, 

6'k = 0, tt/2 ,7r, 37r/2. (3.7) 

In addition, one can easily deduce from Eq. H3.4|l that 

A5 = ^Ask, (3.8) 

k 



^-^Dic^,c)^^^^^D{M ^ COS (^k) CkT + » sin (^k) e-'^'^cli^^ 

^^^Di,p,c)^_^^^^Di4..c) ^ cOS (^k) C^ki " ^ sin (^k) e-^^-'cj^p 

where ^k — I'f'kl f^nd ipy^ — arg((/<k)- Substituting them into Eq. (|3.2(l . we have 

5(0,/?) =ln(Tr(e-^«(^))) +2/3^e(k) [cos\d^)f{e{k))+sm\9u)f{-e{k))] 

k 

I 
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where 



/3£(k)tanhf^4^) [l-cos(2(56'k)] > 0, 6lk = 0, or tt 

^«k=<( )/<^.,{ (3-9) 

-/3e(k)tanhf^j [l-cos(256'k)] < 0, = 7r/2, or 37r/2, 

with (50k being the deviation of 0k f^rom the corresponding solution. 

For the solution that all ^k = or tt, we have AS" > 0, which meets the requirement of Eq. (|2.5|l . thereby, this 
solution corresponds to a stable phase. Observe that the order parameter ^ = when all 0k — 0, one recognizes 
immediately that this phase is just the normal phase of the system, which can also be seen from the following equation. 



13) ^ Tr{F{c)p{H'{cb, c))) = Ty {F (c) p{H (c))) , V0k - 0, or tt. (3.10) 

We are familiar with the properties of this normal phase from the standard books on quantum statistical mechanics. 
What it is worth stressing here is that the normal phase is proved to be stable at any temperature within the framework 
of the extended ensemble theory. Of course, that is reasonable, and just what we expect. 

As to the other solutions, they also have physical senses though they correspond to instable phases because the 
entropies for them are not minimal. Let us first discuss the phase with all 0k = 7''/2 or 37r/2. We find 

T{cb,f3) ^Tr{F{c)p{H'{^,c))) =tAf{c) , V0k = ^, or-, (3.11) 



where p' = l/iksT') with T = -T < 0. That is to say, 
the temperature of the system goes negative. This shows 
that negative temperatures root from the same micro- 
scopic mechanism as that for phase transitions, so we call 
this instable phase the negative-temperature phase. In 
comparison, the normal phase is said to be the positive- 
temperature phase (T > 0). This negative-temperature 
phase can not be realized physically because its internal 
energy is unbounded above. 

Apart from the negative-temperature phase, all the 
other instable phases are partially negative-temperature 
phases, which means that the particles of the system 
are distributed partially among the negative-temperature 
(single-particle) states, and partially among the positive- 
temperature (single-particle) states, 



el3'e('') 



+ 1' 



0k = 0, or TT 

0k 7r/2, or 37r/2. 



(3.12) 



That is to say, some of the states take or tt, the others 
take 7r/2 or 37r/2. Each k with 0k = or tt represents a 
positive-temperature state, and each k with 0k = 7r/2 or 
37r/2 represents a negative-temperature state. Because 
the partially negative-temperature phase is an instable 
phase, it can emit photons when particles transfer from 
a negative-temperature state into a positive-temperature 
one. Therefore, a partially negative-temperature system 
can constitute a laser. 

To make that point more specified, let us consider a 
semiconductor that is described by the following Hamil- 



tonian, 

H{c) 5Z'^"('*)('=",kT'=".kT +cl__k|C„,_ki), (3.13) 



where n denotes the energy-band index, £n(k) the energy 
of the electrons in the nth band relative to the Fermi level 
of the system, and Cn^ko- the annihilation operator of the 
electrons with momentum k and spin <j in the nth band. 
Replacing the D(0, c) of Eq. itT^ with 

D{(t),c) = ^((/'„,kC„,-kiC„,kT + '/>i,kci,kTci,-ki)' 

n,k 

(3.14) 

and following the same procedure as for the ideal Fermi 
gas, one can easily get 



'^LkCT'^n.ko-) — 



^k = 0, or TT 

,k - 7r/2, or 37r/2, 



(3.15) 

where 0„_k = l^n.kl- 

Suppose that there are, for example, two electrons per 
site. If the system is at zero temperature and stays in 
the normal phase where all 0„^k = or tt, the valence 
band, i.e., the 0th band, is fully filled, the conduction 
band, i.e. the 1st band, and all the other higher bands are 
empty. This is the well-known picture for semiconductors 
(or insulators), which is depicted in Fig. 1. 

Now, let us consider such a partially negative- 
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2th band 



conduction band 



valance band 



FIG. 1: Energy bands for the normal (positive-temperature) 
phase of a semiconductor where ^ denotes the chemical po- 
tential. 



temperature phase, 



7r/2, or 37r/2, n = l 
0, or TT, n ^ 1. 



(3.16) 



If this phase is at zero temperature, the conduction band 
is fully filled, all the other bands are empty. It implies 
that the population is inverted betv^feen the valence and 
conduction bands, which is depicted in Fig. 2. As men- 



2th band 



the energies are released by means of emitting photons. 
That is somewhat an ideal semiconductor laser, it is a 
"two-level" system, a more realistic laser can be achieved 
in the so-called "three-level" system, 



0, or TT, 
7r/2, or 37r/2, 
0, or TT, 



n = 

3k for n — 

otherwise, 



(3.17) 



where the 0th and 2nd bands are partially filled, the 1st 
band is empty, which is depicted in Fig. 3. Here, the pop- 
ulation is not inverted between the 0th and 2nd bands, 
i.e., the population of the 2nd band is less than that of 
the 0th band, but it is inverted between the 1st and 2nd 
bands. This kind of "three-level" semiconductor laser has 
been realized in experiments |33| | . For this reason, we call 
the partially negative-temperature phase the laser phase. 
Obviously, a laser phase cannot be described by the orig- 
inal Gibbs ensemble theory 0. 



2th band 



conduction band 



"valance band 



conduction band 



FIG. 3: Energy bands for a laser phase of a semiconduc- 
tor where /i denotes the chemical potential. Here, both the 
valance and 2th bands are partially filled, the population is 
inverted between the 2th and conduction bands. 



valance band 
-R 



FIG. 2: Energy bands for a partially negative-temperature 
phase of a semiconductor where /i denotes the chemical po- 
tential. Here, the population is inverted between the valance 
and conduction bands. 

tioned above, such a population-inversion phase is insta- 
ble, it will finally turn into the stable phase, i.e., the 
normal phase, with the electrons in the conduction band 
hopping into the valence band and simultaneously releas- 
ing energies to the surroundings. It will make a laser if 



Traditionally, negative temperatures are something ob- 
scure, and hard to understand why they arc higher or 
hotter than positive temperatures 0, 0, |36| . From the 
viewpoint of the extended ensemble theory, the reason 
is that both the negative-temperature and laser phases 
are instable, they will trend, by giving up parts of their 
internal energy to the surroundings, towards the stable 
phase, viz., the positive-temperature phase. That is the 
microscopic interpretation about negative temperatures. 

Finally, we note that the ideal Fermi gas cannot pro- 
duce superconductivity because there does not exist any 
Cooper pair in any case of the solutions of Eq. (|3.7(l . 

(c_kiCkT) - Tr(c_kiCkTP(i?'(0, c))) = 0, (3.18) 

this result is just as expected, and right in accordance 
with the BCS theory 0,11 0. 
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In summary, it is proved, within the framework of the 
extended ensemble theory, that the ideal Fermi gas can 
not produce superconductivity, its normal phase is sta- 
ble at any temperature. Besides, the extended ensemble 
theory gives a microscopic interpretation to the negative- 
temperature and laser phases: they originate from the 
same microscopic mechanism as that for phase transi- 
tions, and most importantly they are instable. 



the order parameter, and then to prove that it can sat- 
isfy the requirement of Eq. (|2.5|l . Unfortunately, that 
is rather difficult because it is impossible to calculate 
S{(j), (3) rigorously when H{c) contains an interaction. As 
a result, we have to seek approximations. Before doing 
so, it is worthwhile to give a brief survey to the self- 
consistent mean- field theory due to BCS 0, 0|, which 
can be summarized as follows. 



4. APPLICATION TO BCS 
SUPERCONDUCTIVITY 



4.1. The BCS Mean-Field Theory 



To study BCS superconductivity with the extended en- 
semble theory, one needs first to solve Eq. H2.4|l to obtain 



H{c) =^ Hmf{c) = ^e(k)(4^CkT +c^_kiC-ki) - A ^(c_kiCkT + ^^cl^^) + (4.1) 

k k ^ 

(F(c)) = Tr{F{c)p{H{c))) =^ (F(c)) = Tr {F {c) p {HuAc))) , (4.2) 

where 

A = .g^ (c.kiCkt) = 5^Tr(c_kiCkTP(^^MF(c))). (4.3) 



Eq. 1)4.1(1 and Eq. (|4.3|l constitute a pair of self-consistent 
equations to determine both the mean-field Hamiltonian 
Hmf{c) and the energy gap A. After Hmf{c) is so deter- 
mined, the average of observable F(c) can be evaluated 
according to Eq. (|4.2|) . 

Obviously, if A 7^ 0, the mean-field Hamiltonian 
Hmf{c) is not invariant under the gauge transformation 
of Eq. 1)1. 3|l . i.e., the gauge symmetry is broken with re- 
spect to Hmf{c)- It should be noted that, as shown in 
Eqs. 1)4.2)1 and 1)4.3)) . Hmf{c) must be used in place of 
H{c) when calculating statistical averages. Otherwise, A 
becomes zero, there is no Cooper pair and superconduc- 
tivity, as has been demonstrated in Eq. 1)1.4)) . That is 
the key point of the BCS mean-field theory. This kind of 
mean-field theory is also widely used in studying other 
phase transitions and SSB, Bogoliubov theory of a we akly 
interacting Bose gas being a famous example [stI ISSj . 
Nevertheless, as pointed out by Emch in Sec. 1.1. f of 
Ref. 13], it contains an inescapable paradox. That is, 
the spectrum of the system Hamiltonian Hmf (c) depends 
on temperature, which contradicts the fact that the spec- 
trum of an operator is an invariant property of this oper- 
ator itself and should not depend on temperature. Now, 
for the extended ensemble theory, the normal and su- 
perconducting phases correspond, as stated in Sec. |21 
to the symmetric and asymmetric representations of the 



same system Hamiltonian, respectively, it excludes such 
paradox absolutely. 

Although the approximation of Eqs. ()4.1H4.3)) con- 
tains the paradox, its results are astonishingly in good 
agreement with the experiments. Therefore, we need to 
search a formally identical but essentially different ap- 
proximation within the extended ensemble theory so as 
to interpret the superconductivity anew. 

In language of Green's function (GF), the above ap- 
proximation can be translated into the following formal- 
ism, 



„t \\MF 



Lo + e(k) 



((ckT|c-ki))r 

k 



t^2- [£2(k)+A2]' 

A 

/(u;)Im((ckT|c-ki))r, 



(4.4) 
(4.5) 
(4.6) 



where {{A\B))^^ denotes the retarded Green's func- 
tion defined with respect to 7Jmf(c). Now, Eqs. 
()4.5)) and ()4.6)) constitute a pair of self-consistent equa- 
tions to determine the two complex-valued functions, 
((ckT |c-ki))^^ and A, in contrast to Eqs. (|n|) and lll3|) 
where an operator is involved. Therefore, Eqs. ()4.4H4.6)) 
are formally adaptable to transplanting into the extended 
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ensemble theory. As is well known, BCS superconductiv- 
ity can be explained in terms of these Green's functions 
[J, Q • We shall reestablish them within the framework of 
the extended ensemble theory. Of course, they must be 
redefined, with respect to c) instead of Hmf{c)- 

4.2. Interpretation of BCS Superconductivity 

Now, we devote ourselves to reestablishing a formalism 
similar to Eqs. I|4.4H4.6| . and interpreting the supercon- 
ductivity with the extended ensemble theory. 

Following Eq. H3.2|l . the entropy of a BCS supercon- 
ductor can be written as 

5(<^,/3) =ln(Tr(e-''^(^))) 

/3^£(k)(4^dk| + d^_ki'^-ki) 
k 

-/35EE4'T^-k'i'^-ki'ikT, (4.7) 

k k' 



where 

dk| = e-*-°('^'=)ckTe''°("*''=) (4.8a) 

(i_ki - e-'-°('^^")c_kie*^('^''=), (4.8b) 

and 



A{c)^Tr{A{c)p{H{c))). (4.9) 

As the usual mean-field approximation, we decouple the 
last term on the right-hand side of Eq. (|4.7|) as follows. 



-k'i ^ f^-kiC^kT = C^k'T -k'i ^ '^-kirfkT • (4.10) 
It results in 



S{(b,(3) = ln(Tr(e-''^(^))j +/3^e(k) [cos2(0k)(4^CkT + cl^iC-ki) + sin2(0k)(ckT4T + c-kic^_ki; 

k 

- i^5E^i^(20k')e^'''''(cL^Ck'T -c-k'iclk'i) x ^ sin(2^k)e-'^''(4^CkT " c_kiclk|), 

k' k 

where Eqs. H4.8|l . (|3.3|l . and p.4|) have been used. 

Substitution of Eq. H4.11|l into Eq. (|2.4|) gives the equations of order parameter. 



2e{k) sin(26'k) + g cos(26'k) ^ '^^^('^k - v^k') sin(26'k')(Ck'T'^k'T - c_k'iC^k'i) = 



sin(26'k) E ^i^('^k - (ySkO sin(26'k')(ck'T'^k'T - c^u'icL^'i) = 0. 

k' 

I 



(4.11) 



(4.12) 
(4.13) 



They have a trivial solution. 



= 0, 



and a nontrivial solution, 



^ fa ^ constant 
e(k) sin(20k) - Acos(26'k) = 0, 



where 



A = - 



kiCkt/ 



(4.14) 



(4.15) 



(4.16) 



The trivial solution always remains zero, it does not 
change with temperature. In contrast, the nontrivial one 
will depend on temperature, it may be zero for some tem- 
peratures, and nonzero for the others. Without loss of 



generality, we shall take ipo — 7r/2, the nontrivial solution 
is then simplified as 



e(k) sin(26'k) - Acos(26'k) 
A = .gE (c-kiCkt) , 

k 



0, 



(4.17) 
(4.18) 



they are a pair of coupled equations to determine both 
6'k and A. 

If g = 0, Eq. l|TTH|l shows that A = 0, then Eq. (|CT7|) 
comes back to the equation (|3.6|l for the ideal Fermi gas. 
Physically, the BCS Hamiltonian of Eq. p. II) can de- 
scribe only the weak-coupling superconductivity. For the 
weak-coupling superconductivity, the coupling strength 
g is quite small, from Eqs. H4.17|l and H4.18|l it follows 
that A and 9u are also small. This makes it easy for us 
to do further approximations, as will be seen from below. 
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In order to solve Eqs. (|4.17() and H4.18|l . one needs where {A,B} denotes the anticommutator of A and B. 
(c_k^Ck|), it can be obtained by the corresponding re- Expanding H'{(j),c) into the power series of 0, 
tarded Green's function ((ckj |c_ki))w defined with re- 
spect to H'{(j},c) ^3^. As is well known, ((ckj |c_ki))[^ 
satisfies the equation of motion, 

w((ckT|c_ki>>^ = ({ckT,c_ki}) 

-((ckT|[c-ki,i?'(0,c)]))^, (4.19) 



c) = H{c) + z c),H{c)] + - [D{cl,, c), [D{^, c),H{c)]] + ■■■, 

and substituting it into Eq. 14.19|1 . we have 

w((ckT|c-ki))c^ = ({ckT,c_ki}) - ((cktl [c-ki,H{c)]))^ - «((ckTl [c-ki, [D{(l),c), H{c)]]))^ 



(4.20) 



2! 



((ckT|[c-ki,[i?('^,c),[i^(0,c),i/(c)] 



(4.21) 



As stated above, is small for the weak-coupling superconductivity, it is thus rational to maintain only the zeroth- 
order term, 



w((ckT|c_ki))^ = -e(k)((ckT|c-ki))^ +5((ckT|CkTX!'^-k'iCk'T)). 

k' 

I 



(4.22) 



With regard to the second GF on the right-hand side, we 
do the factorization as in Eq. (|4.1Q|1 . 

9^{c--k'icu'^) = g^d-u'idk'^ = A, (4.23) 

k' k' 

this leads to 

[a; + e{k)] ((ckt |c-ki». - A((ckT (4.24) 

The GF on the right-hand side can be obtained by the 
same procedure as for ((ck| |c_kj.))i^, 

[u; - £(k)] ((ckTl4^»c. = 1 + At((ckT |c-ki))c., (4.25) 
where 



From Eqs. I4.24H4.26|) . it follows that 

((cktlckt))" = 
((ckT|c-ki))t^ = 



t XX ^ ^ + g(k) 

T''" t^2_ [^2(k)+A2]' 

A 



A = -<?E / — /(w)Im((ckT|c-ki)). 

1. ^ — OO ^ 



(4.26) 

(4.27) 
(4.28) 
(4.29) 



They form a self-consistent mean-field solution to the 
Green's functions responsible for BCS superconductivity 



within the framework of the extended ensemble theory. 
Comparing them with Eqs. (|4.4H4.t)|) . one sees that the 
present solution is identical to the BCS solution except 
that the statistical average is now defined with respect 
to iJ'((/), c). Thus far, the BCS mean-field results have 
been transplanted into the extended ensemble theory, the 
paradox removed. 

With the help of Eq. W^ . Eqs. (ji:^ and (|CT7|) 
can be simplified as 



A = 



7k = 



1 tanh (^/3Ci 

2^ ^ 



iarcsin(A) , 
-iarcsinfA 



e(k) > 
e(k) < 0, 



where 



^k = v/e'(k)+A2. 



(4.30) 
(4.31) 

(4.32) 



Eq. H4.31|l together with Eq. H4.30|l gives the mean-field 
solution to (0k = i^k), the order parameter for BCS 
superconductivity. 

_Eci. (|4.30() is familiar in the BCS mean- field theory 
ic nontrivial solution, 

A = [ ^«™' lll^ (4.33) 
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where is determined by the equation, 

_ 1 ^ tanh (|/3c£(k)) 
'"2^4^ £(k) 



(4.34) 



with 13c = l/ksTc. From Eqs. and l|OT|l . we find 



zero, T > Tc 
nonzero, T < Tc- 



(4.35) 



Since (j) is the order parameter of the system, Eq. (|4.35|l 
indicates that Tc is the transition temperature for BCS 
superconductivity. 

Compared with Eq. I|4.35|l . Eq. (|4.33|) shows that the 
energy gap A behaviors with temperature just hke an or- 
der parameter. However, this behavior of the energy gap 
is not a property of itself, but derives from the property 
of the order parameter, as can be seen from the equation. 



A = .9 ^ Tr(c_uxCkTP {H'{cf,, c)) ) = | 



zero, (f> ~ 
nonzero, 7^ 0, 

(4.36) 

which demonstrates that the forming of Cooper pairs is 
a direct consequence of the gauge-symmetry breaking 



caused by the internal spontaneous field (f>. In the ex- 
tended ensemble theory, the order parameter is a physi- 
cal quantity more fundamental than the energy gap. This 
picture is significantly distinct from the BCS mean-field 
theory where the energy gap is the most fundamental 
quantity that is responsible for the superconductivity. 

In sum, Eqs. g^B, EiSHIl, 6301), and l|OT|l consti- 
tute a self-consistent mean-field theory for BCS supercon- 
ductivity within the framework of the extended ensemble 
theory, its nontrivial solution reproduces the BCS mean- 
field results on the conventional low- Tc superconductiv- 
ity. 

We are now confronted with the important task to 
prove the stability of this nontrivial solution. By expand- 
ing the S{<p, P) of Eq. H4.11|l into a Taylor or Volterra 
series in the neighborhood of the nontrivial solution 0o, 
one finds 

AS ^ S{(l),l3)-S{(j)o,P)^d'^S+d^S+5'^S+--- , (4.37) 



where (5" 5 represents the nth power term of the expan- 
sion of S{(j),f3), or the nth variation of S. Specifically, 
the second variation 6'^S has the form, 



x2q 1 a ,2 tanh(l/3ek) tanh(l/3a') . \^^^rX^c f w ^ At \txa ^2 

S S ^ -f3gA 2^2^ 7^^ [Sifik - Sifii,,) +2/32^Ckitanh(-/3^ki)(56'ki) 



k k' 



tanh(-/3ekj [Se^.f - E 



tanh(i/3aj ,w^tanh(i/3a^) 



k2 k; 



Ck 



6 



(4.38) 



r 



where Sipk and SOk represent the variations of (fk and 9^ 
from the nontrivial solution c/jq, respectively. Here, the 
set of k has been separated into two subsets: the subset of 
ki and the subset of k2, where ki and k2 satisfy e(ki) = 
and e(k2) ^ 0, respectively. All the contributions from 
the subset of k2 are included by the term within the curly 
brackets. 

Let us prove that this term is definitely positive, it can 
be done by proving that the eigenvalues of the matrix M, 



k.k;, =2fk. tanh(-/3aj'5^,k; 



tanh(i/3a2)tanh(i/3Ck' 



Cka Ck^j 

are all positive. The eigenvalue function for M is 



\n{X)\ = o, 



where 



f^k.k;(A) 



k2k. 



M, 



kgk. 



(4.39) 



(4.40) 



(4.41) 



Observe that the second term on the right-hand side of 
Eq. (|4.39|) is a product of the two factors corresponding 
to the row k2 and the column k2 respectively, we can 
reduce Eq. 14.4011 into 



where 



^^k.k;(A) 



Obviously, 



MX)\ = o, 



A -2a, tanh(i/3aj 
g [£(k2)tanh(i/3aj/a2 



2'^k2k; 



A = 2a. tanh(-/3aj >0 



(4.42) 



-1. (4.43) 



(4.44) 



is an eigenvalue of Eq. (|4.42|) because at least two rows 
(±k2) of the determinant |ri(A)| are identical. Apart 
from those eigenvalues, there are possibly other ones. 
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they satisfy the equation, 



[£(k2) tanh(l/3CkJ/ek.] 
2^k, tanh(i/3aj- A 



(4.45) 



To derive this equation from Eq. H4.42I) , it is sufficient to 
heed that every element of f2(A) contains the number 1. 
If A < 0, one has 



[£(k2)tanh(i/3aj/ek.]' 1 ^tanh(i/3a) 



< 



tanh(i/3ekj - A 2 



Eqs. H4.45|l and (|4.46(l lead us to 



(4.46) 



This inequality contradicts the result of Eq. 14.3011 , 



1 ^ tanh(i^a) 
.9 2^ ^ < 1- 



2 V 



(4.48) 



That is to say, the root of Eq. 14.45() must be greater 
than zero if it has any. Combination of this result with 



Eq. H4.44|l demonstrates that all the eigenvalues of the 
matrix M are definitely positive, which ends our proof 
for the positivity of the term within the curly brackets of 
Eq. Km . 

This positivity implies that 



S^S > 0. 



(4.49) 



When S'^S > 0, we know from Eq. l|137jl that AS" > 0, 
it meets the requirement of Eq. H2.5|l . However, if 5^S = 
0, we must generally examine 5^S, 6^S, or even more 
higher variations of S to check whether AS' > 0, which 
is obviously complicated to handle. Fortunately, we need 
not do that in the case considered now. When 6^S = 0, 
Eq. shows that 



S9u 



0, if T > Tc, and 



(56'ki =0 if T < Tc. 

sei,, = 0, 

(4.50) 

With them, one can easily verify that S{(j), 13) — S'((/)o, /3). 
Namely, AS* = if = 0. 

To sum up, AS* > iiS^S > 0, and AS* = iiS^S = 0. 
Since S^S > 0, we conclude that the requirement AS" > 
is satisfied, the nontrivial solution is stable. 



As regards the trivial solution, (/> = (A = 0), it is identical to the nontrivial one when T >Tc, and thus is stable 
at T > Tc. If T < Tc, we find 

6^S = fi ||]2e(k2)tanh Q/3e(k2)) (50kj' - 5 1] tanh Q/3£(k2)) tanh ^/^^(k'a)) ^^k.-^^k^j ■ (4.51) 



Accordingly, it has positive eigenvalues. 



1 



A = 2£(k2) tanh ( -(ie{\^2) ] > 0, (4.52) 



In addition, we also have 



[tanh(l/3£(k2))]^ 
2£(k2)tanh(i/3e(k2)) - A' 



(4.53) 



In contrast to Eq. H4.45|l . this equation has a negative 
root of A when T < Tc. To see it clearly, let us consider, 
for instance, the zero temperature case, 

f*hujr 



1 Pllulo 1 



A' 



e(k)>0 

(4.54) 

where, as usual, A/'(0) denotes the density of states at the 
Fermi level, and ujd the Debye frequency. Apparently, it 
has the solution, 



A 



2nw 



D 



exp(^Ar(oy) - 1 



< 0. 



(4.55) 



The existence of both the positive and negative eigenval- 
ues implies that the trivial solution is a saddle point at 
T < Tc. Therefore, the term within the curly brackets of 
Eq. H4.51|l is indefinite, that is, it can change in sign. As 
a consequence, the trivial solution, or rather the normal 
phase, is unstable at T < Tc. That is just the Cooper 
instability HHQ. 

Of the two solutions, the nontrivial is the only one that 
is stable within the whole temperature range. 

With the nontrivial solution obtained, other physi- 
cal quantities can be easily expressed and calculated 
in terms of the two Green's functions ((ckt |c]^t))'^ ^'^d 
((ckt|c_ki))w, as was done in Refs. 0, and 0|, 
with the results unchanged and the same as the BCS 
theory. 

So far, BCS superconductivity has been interpreted 
within the extended ensemble theory. In this interpre- 
tation, the order parameter for BCS superconductivity 
evolves with temperature according to the principle of 
least entropy. At high temperatures (T > Tc), it is zero 
and stable, the system Hamiltonian realizes the represen- 
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tation with perfect gauge symmetry, the resulting phase 
is normal and disordered. As temperature decreases and 
goes below the critical temperature (T < Tc), the zero 
solution becomes unstable, instead, there arises a new 
stable solution of the order parameter, it is nonzero. At 
the same time, the system Hamiltonian transforms into 
a new representation with gauge symmetry broken, the 
electrons are hence formed into Cooper pairs, and the 
resulting phase gets superconducting. Thus and so, the 
extended ensemble theory answers why, when and how 
the electrons can coagulate to form Cooper pairs and su- 
perconductivity, the problem posed by Born and Fucks 

Lastly, it should be noted that the BCS and Landau 
theories are unified into a single formalism within the 
framework of the extended ensemble theory. 



the hypotheses in Sec. 12 to consider the phase-transition 
operator, 

D{v,b)=J2ivlbi. + Vi.bl), (5.5) 

k 

where the internal field rj is the order parameter for BEC. 
The reason for introducing Eq. (|5.5|l can be seen from 
comparing Eqs. (|5.4|l and H5.5|l with Eqs. I|1.4(l and H2.2() 
respectively. 

We shall consider first the finite volume case where V 
and TV are both finite, and then the thermodynamic limit 
case where V +oo, and N +oo, but N/V is finite. 
As will be seen, the two cases are significantly different. 

5.1.1. The Finite Volume Case 



5. BOSE-EINSTEIN CONDENSATION 

After the discussion of BCS superconductivity, we pro- 
ceed now to study Bose- Einstein condensation (BEC). 
We shall concentrate on three systems: the ideal Bose 
gas, the photon gas in a black body, and the ideal phonon 
gas in a solid body. They are the most important cases 
of Bose systems: the first is a system with the conserva- 
tion of particles, but the other two are not; the second 
belongs to gauge fields whereas the third does not. In 
addition, we shall also discuss the quantization of Dirac 
field from a standpoint of the extended ensemble theory. 

5.1. The Ideal Bose Gas 

Let us begin with the ideal Bose gas. Its Hamiltonian 
reads as follows, 

ff(&) - - M)6lc&k, (5.1) 

k 

where 6k (^k) is the annihilation (creation) operator for 
the bosons with a momentum fi,k, and ek = /i^k^/ (2m) 
the single-particle energy, and ji the chemical potential. 
This Hamiltonian is gauge invariant, 

G{d, h)H{b)G\d, h) = H{b), (5.2) 

where 

G{^,b) = e-'^^'^''^''^, d e [0,27r). (5.3) 

The invariance has an important consequence, 

(6k) = Tr(6kp(i?(6))) = 0. (5.4) 

It signifies that the condensation amplitude will be zero 
forever if the gauge symmetry does not break down. 

In order to examine whether the gauge symmetry can 
break down spontaneously or not, one needs, following 



In this case, the wave vector k is discrete. The entropy 
of the system can be written as 

5(77, 13) = S{0, /?) + /? ^(ek - (5-6) 

k 

where 

5(0, /3) = -Tr( In (p(i/(6))) p{H{b))) (5.7) 

is independent of the order parameter ry. 
According to Eqs. (|2.4|l and H2.5(l . we have 

1^ = 0, (5.8) 
A5 = /3 ^(ek - t^)Sri5r^k > 0, (5.9) 

k 

where (Jyyk represents the variation of rjk from the solution 
given by Eq. (|5.8I) . Eq. (|5.9|) shows that the chemical 
potential must satisfy the following condition, 

M < ek. (5.10) 

As Ek > 0, this condition is equivalent to 

/i<0. (5.11) 

This inequality is the physical condition for the stability 
of the system. 

In combination with Eq. H5.6|) . Eq. (|5.8|l gives us the 
equation of order parameter, 

(ek-M)'7k = 0. (5.12) 

It has two solutions: the trivial one, 

r?k = 0, (5.13) 

and the nontrivial one, 

= '7o<5k,o, only if ^ = 0. (5-14) 
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For the nontrivial solution, /i = 0, it causes a contra- 
diction, 

= 7V-^Tr(5[5kp(i?(5))) 

k 

= ^-j:^p^^ (5-15) 

k 

where the left-hand side can not be negative, but the 
right-hand side is negative because the term of k = 
is infinite. Therefore, the nontrivial solution must be 
discarded. 

With regard to the trivial solution, it corresponds to 
the normal phase of the system. It does not cause any 
contradiction, 

N ^ ^Tr(b[b,p(g(b))) ^ J2 mJ,, , r (^-'^^ 
k k 

this equation can be satisfied at any temperature. In 
fact, it is just the equation for the chemical potential /i, 
its solution will give = Because N is finite, fi 

will be less than zero irrespective of temperature, i.e.. 



where 



fi{T) < 0, for T > 0. 



(5.17) 



It implies that the normal phase of the system is always 
stable, 



AS > 0. 



(5.18) 



This proves that there is no phase transition in the finite 
volume case, the system will stay in its normal phase 
forever. 



5.1.2. The Thermodynamic Limit Case 

In this case, V +oo, we should use the density of 
entropy instead of the entropy itself. 



slri,j3) = lim 



V^+oo V 

s(0, 13) + /? hm i ^(ek - (5-19) 



5(0, /?) = - hm -Tr( In (p(i/(6))) p{H{h))) 

V—*+oo V 



(5.20) 



is a term irrelevant to rj. The limit V +oo can be 
achieved as follows. 

As usual, suppose that the system is box normalized. 



k 



277 



x{y,z) 



x{y,z) 



(5.21) 



where Z denotes the set of integers, and Ix, ly, and Iz 
the dimensions of the box along the x, y and z direc- 
tions respectively {V = Ixlylz)- As is well known, each 
k corresponds to an eigenfunction, all the eigenfunctions 
constitute a Hilbert space. Here, the Hilbert space is 
a complete space of integrable functions, with the in- 
ner product defined as a Lebesgue integral. It should 
be pointed out that the inner product must be defined 
as a Lebesgue rather than Riemann integral because the 
former can ensure the completeness of an integrable func- 
tion space whereas the latter can not .42\. This implies 
that all the problems relevant to a Hilbert space of inte- 
grable functions must be treated with Lebesgue theory. 
On the other hand, Lebesgue measure and integration 
is itself the mathematical foundation of probability the- 
ory and statistical physics . We shall there- 
fore perform the thermodynamic limit according to the 
Lebesgue theory of integration. 

Formally, Lebesgue integration includes both the sum 
over discrete numbers and the integral on a continuous 
region. When V is finite, a sum over k is a Lebesgue in- 
tegral of discrete form. It will transform into a Lebesgue 
integral of continuous form as V ^ +oo. To show the 
transformation, let us consider, e.g., the sum over k in 
Eq. (|5.19(l . In the first place, we write it into an explicit 
form of Lebesgue integral. 



^ - (,^1') k 



/i(p)dp. 



(5.22) 



where 

m3 



(— (X),+oo) X (— oo,+oo) X (— cx3,-|-oo), (5.23) 



m(ylknM^) 



V 



(5.24) 



27r 2tt 271 

Aw = [kx,kx + ^) X [ky,ky+ — ) X [kz,kz + — ), (5.25) 



(5.26) 



As usual, XAk(p) represents the characteristic function 
of the measurable set Ak, and m(Ak H M^) the Lebesgue 
measure of the set Ak- From Eqs. (|5.25H5.27|l . it follows 
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that Now, in the second place, let us take the thermodynamic 

limit on both the sides of Eq. I|5.22|) , 

lim h{p) = (cp - /i)??],?7p 

V — »+oo 

= [e{p)-^lWiPMp)■ (5.28) 



lim y^fsk - /-i)?7l?7k = — ^ lim / /i(p)dp = — ^-^ / lim /i(p)dp = — ^-^ / [e(k) - d77^(k)77(k)dk. 



r 



(5.29) 



In such a natural way, a discrete sum over k will trans- 
form into a continuous Lebesgue integral in the thermo- 
dynamic limit. 

With the help of the above equation, one arrives at 

[e(k) - ^^] r^ik) = 0, (5.30) 

As = /3-i-3 / [6(k) - iASvH^)Sv{k)dk > 0. (5.31) 
(/tt) J 

The second equation means that the chemical potential 
must be less than or equal to zero. 



/i<0. 



(5.32) 



it is the condition for the stability of the system. The 
first equation is equivalent to 



[e(k)-M],7t(k)r;(k)=0, 



(5.33) 



which can be easily verified with respect to the two cases: 
?7'f(k) = and r]'^ (k) 7^ 0. Eq. has two solutions: 

the trivial one. 



V{k) = 0, 



and the nontrivial one. 



77(k) = eVW, only if A* - 0, 



(5.34) 



(5.35) 



where ^ is a complex which does not depend on k, and 
i5(k) the Dirac 5 function. 

For the trivial solution, the system is in its normal 
phase, this phase can exist only at high temperatures, 
which can be easily deduced from the equation of chem- 
ical potential, 

n = lim — 

V'^+co V 



V — >+oo V 



1 



(27r)^ J e'5(^(k)-p) _ I 



dk 



0+ 



p/5(e-M) — 1 ' 



(5.36) 



where we have reduced the Lebesgue integral into an im- 
proper Riemann integral in the end. Because, as shown 
by Eq. (|5.32() . /i must be less than or equal to zero, 
the above equation can not hold if T < where is 
determined by 



_ (2m)3/2 1-+°- eV^ 



As is well known, this gives 

2TTh^ ( n 



mkB VC(3/2 



2/3 



(5.37) 



(5.38) 



where C(x) denotes the Riemann zeta function. 

For the nontrivial solution, /x must be zero, the order 
parameter is determined by 



1 



(2^)^ 



?7^(k)r/(k)dk^ 



1 



1 



{2T:f J e/3^(k)-l 



dk 



1 . (2m)3^ .1/2 



(2^)^"^^' 47r2n3 



de- 



(5.39) 



where we have used the fact that the contribution from 
the state of k = is zero when reducing the Lebesgue 
integral into the improper Riemann integral in the end. 
Because > 0, the above equation can not hold if 
T>Tc. 

From those discussions, it follows that 



zero, T > Tc 
nonzero, T < T„ 



(5.40) 



That is to say, a phase transition will happen at T = Tf,. 
Because As > when /i < 0, the normal phase (77 = 0) 
exists and is stable a.t T >Tc, the ordered phase (77 7^ 0) 
exists and is stable at T < Tc- Here, we note that it is 
a natural consequence of the extended ensemble theory 
that the chemical potential of the ideal Bose gas is fixed 
at zero = 0) in the ordered phase (T < Tc), in contrast 
to the Einstein's treatment where it is fixed by physical 
insight [illiell. 

Obviously, the above results are the same as the BEG 
given by Einstein if ^^^/(27r)^ were regarded as 
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the density of the bosons condensed onto the state of 
k = 0. Therefore, the transition happening at T = 
must be identified physically as the Bose-Einstein con- 
densation. It should be pointed out that, generally, the 
quantity can not be interpreted as the density 

of condensed bosons, as will be clarified in Sec. 17.11 

So far, we have presented a new description to the 
Bose-Einstein condensation of the ideal Bose gas. In 
comparison with the Einstein's description |45l |46j , we 
now describe the BEC with two new conceptions: spon- 
taneous symmetry breaking and Lebesgue integration. 
Within this new description, the BEC happens simul- 
taneously with the spontaneous breaking of the gauge 
symmetry; the integration is defined and performed in 
Lebesgue way. That is fundamentally different from Ein- 
stein's description ji^ where the BEC happens with 
no symmetry breaking, and the integration is done in 
Riemann way. 

Here, it is worth emphasizing the important role played 
by Lebesgue integration in the study of the BEC of the 



ideal Bose gas. As has been seen, there arises the Bose 
distribution function, 

= (5.41) 

which is unbounded above on the set E. Because the 
integrand of a Riemann integral must be bounded ev- 
erywhere, such a function is not Riemann integrable and 
can not be handled with Riemann integration. Fortu- 
nately, that makes no trouble. As pointed above, it is, in 
fact, unnecessary for us to handle the Bose distribution 
function with Riemann integration; what really confronts 
us in statistical physics is Lebesgue other than Riemann 
integration. For Lebesgue integration, it permits more 
general functions as integrands, and treats bounded and 
unbounded functions on an equal footing. In the sense 
of Lebesgue integration, the Bose distribution function is 
integrable, which can be shown as follows. 



\ ^ ]^ Ifl l/'l l/'l 

y^™oo V ^ e^^^k _ 1 " (27r)^ J ef^^^) - l^'^ ~(27r)^ J et^^^) - l'^''^(27r)^ J e'^'C') - l'^'^ 

dk^iS^/^ d.-^, (5.42) 



(27r)^ J eP<^)-l An'^h^ 



where the set i?o = {k = 0} includes one point only. Eq 
is a null set, i.e., its measure is zero, thus, the integral 
on it vanishes. The rest part, i.e., the integral on R'^\ii^oi 
reduces finally into an improper Riemann integral. As 
is well known, the integral on a null set always vanishes 
irrespective of the integrand values on this set, even if 
they are infinite. Thereby, the integral on any single 
state equals zero, the ground state (k = 0) is not more 
advantageous than the other states (k 7^ 0). In other 
words, the singularity of /(e (k)) on the point k = has 
no special meaning, both in mathematics and in physics. 

This aspect makes the present treatment quite distinct 
from that of Einstein, where the integration is handled 
in Riemann way, with a special treatment given to the 
ground state of k = 0. Intuitively, Einstein thought that 
the bosons would gather in the ground state if T < T^ , 
and thus the contribution from this state is nonzero |45l 
li^ . that is, 

n ~ lim — —3 

k 

f2rnl3/2 /•+00 1/2 

where no > is the particle density gathering in the 
ground state. Comparing Eqs. (|5.39|) . (|5.42|) and (|5.43l) 



shows more clearly the difference between the present 
treatment and Einstein's. From this difference it follows 
that the quantity ^'''^/(27r)'^ needs not to be interpreted 
as the density of the bosons condensed onto the state of 
k = 0. In fact, Eq. (|5.39|l is physically just an equation 
of order parameter. Here, it should be noted that the 
Bogoliubov theory of weakly imperfect Bose gas [s^, 113 
is based on Eq. (|5.43|) rather than Eq. H5.42|l . At last, 
why did Einstein not use Lebesgue integration? there 
is a historical reason: Until 1925 when Einstein |45L |4^ 
studied the BEC, Lebesgue integration had not yet been 
introduced to probability theory by Kolmogorov, who did 
it in 1933 0. 

Also, it is significant to compare the two Lebesgue in- 
tegrals represented in Eq. (|5.15|) and Eq. (|5.39|l : Both 
integrands are unbounded on R"^; however, the former is 
infinite on a set with a nonzero measure (m(Ak=onM'^) = 
(27r)'^ IV\ the latter is finite almost everywhere, or infi- 
nite just on a null set (to {Eq n M'^) = 0). This difference 
leads to the result that the nontrivial solution must be 
discarded in the finite volume case whereas it has phys- 
ical sense in the thermodynamic limit. In a word, the 
ideal Bose gas can produce BEC in the thermodynamic 
limit, but it can not in the finite volume case. That is a 
rigorous conclusion deduced from the extended ensemble 
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theory. 

As another instance, the ideal Bose gas shows unam- 
biguously that Gibbs ensemble theory needs some kind of 
extension. Otherwise, the conservation of particles will 
be broken when T < Tc, as can be seen by comparing 
Eqs. ifO^ and (|0^ . 

As the only one of exactly solvable quantum phase 
transition, the BEC in the ideal Bose gas demonstrates 
that the extended ensemble theory postulated in Sec. [3 
is feasible in physics for us to describe phase transitions. 

Lastly, we would like to stress that the condition of 
Eq. H5.32|l . i.e., /i < 0, is not a mathematical limit but 
just a physical requirement. Like the case of < 0, the 
extended ensemble theory is also well defined mathemat- 
ically for /i > 0, as is shown in Appendix Remov- 
ing the limitation on the chemical potential is physically 
reasonable, it enables the conservation of particles to be 
handled freely, and makes it easy to describe the BEC in 
the interacting Bose gas, as can be seen in Sec. 17.11 

Remark 1: The ideal Bose gas is a touchstone for the 
theory of phase transitions because it is the only one of 
the many-body systems that can show a quantum phase 
transition and can be solved rigorously. Any theory of 
phase transitions must be examined by the ideal Bose 
gas. It can not be accepted in physics if it can not give 
a rigorous solution to the ideal Bose gas; and it can be 
further applied to the interacting Bose gas only after it 
has given a rigorous solution to the ideal Bose gas. 

Remark 2: In 1926, Uhlenbeck raised an objec- 
tion to Einstein's viewpoint on the condensation, he ar- 
gued that the ideal Bose gas can not produce BEC if 
the volume of the system remains finite. This objection 
caused an exciting debate during the van der Waals Cen- 
tury Conference in November 1937 Evidently, the 
present theory has given the exact answer to the debate: 
Uhlenbeck and Einstein are both right, which is also in 
accordance with the agreement reached finally in the con- 
ference. 

Remark 3: In the momentum space used here, the 
condensed phase is shown to be homogeneous. In Sec. 
17.11 using the real space, we shall prove further that the 
condensed phase can only be homogeneous. There can 
not exist any supercurrent or quantized vortex in the 
ideal Bose gas. 

Remark 4: The rigorous solution for the ideal Bose 
gas shows that the extended ensemble theory holds in the 
critical region. 

5.2. The Photon Gas in a Black Body 

For the photon gas in a black body, its Hamiltonian 

reads, 

H{a) = 51 ( aL^k. + ^ ) , (5.44) 

where Coulomb gauge is used, Okjy (^kjy) destruc- 
tion (construction) operator for the photons with mo- 



mentum k and polarization u, and Wk = ^ck represents 
the photon energy with c the speed of light. 

Analogous to the ideal Bose gas, this system has the 
phase-transition operator, 

m a) = E(^k-°k. + Ck.aL), (5-45) 

ki/ 

where ^ is the order parameter for BEC. 
From H{a) and D{£^, a), it follows that 

wk^k. = 0, (5.46) 

AS ^ PY.'^^^^U^^'- (5-47) 

k,cr 

The first is the equation of motion of order parameter, it 
has the solution, 

^k,. = KSk,o, (5.48) 

where Xi, is a complex constant which does not depends 
on T. Owing to the fact: Wk > 0, this solution is stable, 

AS" > 0. (5.49) 

To determine Aj/, it is sufficient to consider the initial 
state at /? = (T — > +oo). Physically, any system must 
stay in its disordered phase at /3 = 0. Therefore, must 
be zero at /3 = 0. Since Xi, does not depends on /3, it is 
always zero. That is to say, 

^ (/?) = 0, for /? > 0. (5.50) 

This equation shows that the photon gas in a black body 
can not produce BEC, and must always stay in its normal 
phase. Eq. (|5.49|l demonstrates that the normal phase 
is stable. This result is also valid for the thermodynamic 
limit case. 

In short, there can not arise BEC in the photon gas in 
a black body, the system stays in its normal phase for- 
ever, Planck's radiation law holds at any temperature. 
Obviously, this conclusion agrees completely with the ex- 
periment. 

5.3. The Ideal Phonon Gas in a Solid Body 

As the third system, let us consider the ideal phonon 
gas in a solid body, whose Hamiltonian reads, 

H{d) = J2 '^k. (dlj^, + , (5.51) 

where dku {dl^^r) destruction (construction) opera- 

tor for the phonons with momentum k and polarization 
cr, WkCT = ^.c^k represents the phonon energy with the 
velocity of sound. 

By comparing Eq. H5.51|l with Eq. H5.44|l . one con- 
cludes at once that the ideal phonon gas in a solid body 
can not produce BEC either, and must also stay in its 
normal phase forever. 
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5.4. Quantization of Dirac Field 

At last, let us consider the quantization of Dirac field 
from the standpoint of the extended ensemble theory hy- 
pothesized in Sec. [3 If one quantizes Dirac field using 
Bose-Einstein statistics, he has the following Hamiltonian 

H{b- d)^J2''p i^lJp^^ - dladp..) , (5.52) 

p,cr 

where p and a denote momentum and spin, respectively; 
bp cr (6p a-) dpcr (dpa) ^-^^ bosonic destruction (con- 
struction) operators; and 

cj(p) = ^J{cpf + {mc^f (5.53) 

with m and c being the electron mass and the speed of 
light, respectively. 

Evidently, this system has the phase-transition opera- 
tor, 

D{v, b; d) = ^{ill^hp^ + S,l„dp„ 

p.cr 

+ Vpahl,+£.pad\,,), (5.54) 

where "q and ^ are the order parameters for BEC. 
Eqs. (|2.4|) and (|2.5|) imply that 

ripa = 0, (5.55) 
^p. = 0, (5.56) 
A5 /3 ^ Up{5r,lJr^p^, - S^U^p^,). (5.57) 

p.cr 

Eqs. (|5.55() and (|5.56() show that the system should al- 
ways stay in its normal phase, but Eq. (|5.57|) demon- 
strates that this normal phase corresponds to a saddle 
point, and is thus instable. Therefore, if Dirac field 
were quantized using Bose-Einstein statistics, the sys- 
tem would be instable and collapse. In other words, it is 
not admissible to quantize Dirac field using Bose-Einstein 
statistics. This conclusion is in accordance with the spin- 
statistics theorem in quantum field theory jl^ • 

On the other hand, if Dirac field is quantized using 
Fermi-Dirac statistics, its Hamiltonian will read as fol- 
lows 

H{h- d)^J2 ^pi^l^bp,, + dl^dp^a), (5.58) 

p,tT 

now bp rj (&p cr) ^'^d dp„ (dp^) are fermionic destruction 
(construction) operators for electrons and positrons. Ac- 
cording to Sec. 121 its phase-transition operator is 

D(0, b; C, d) = ^[(0kfo-ki6kT + ^l&lcT^-ki) 

k 

+ (Ckd-kirfkT + CUi^dl^^)], (5.59) 

where cj) and ^ are order parameters. Following the ar- 
guments of Sec. 131 one recognizes immediately that the 



system is always stable in its normal phase. This holds 
even in the original representation of the Hamiltonian 

p s=3 / 

where 

( bp^s, s = l,2 

Dicj), C, C) = ^[(?!'kC_k,2Ck,l + '/'k4,icLk,2) 
k 

+ (CkC-k,4Ck.3 + CIcIs^Ua)] ■ (5.62) 

In contrast to the electron-positron representation of Eq. 
(|5.58|) , a minus sign appears in the original representation 
of Eq. (|5.t)U|) , but it can not infiuence the stability of the 
system. Anyway, there is no problem when Dirac field is 
quantized using Fermi-Dirac statistics. 

In this section, the systems concerned are free boson 
fields, they can all be solved rigorously. We shall proceed 
to the interacting systems in Sec. [T] Generally, they can 
not be solved rigorously, approximations are needed. 



6. PHYSICAL ORIGINATION FOR PHASE 
TRANSITIONS 

Before proceeding forward to the interacting Bose sys- 
tems, it is worth making a research into the physical orig- 
ination for phase transitions. 

Again, let us take BCS superconductivity as an in- 
stance. According to Eq. H2.3|l . the statistical average of 
a one-body operator can be reformulated as 

J dr(^t(r)/(r)^(r)) = J drTr(^t(r)/(r)^(r)p(iJ(c))) 
= J drW{mr)W), (6.1) 

where ip{r) is the electron field, /(r) stands for the one- 
body operator in Schrodinger picture, and 

V^(r) = e-*-°('^''=V(r)e'-°^'^''=^ (6.2) 

After the phase transition {(j> ^ 0), the electron field ip{r) 
is separated into two fields, 

V^(r) = V(r) +¥'(!■), (6.3) 

where 

ip{r) = e~'-°("*''=V(i-)e'-°^'^''^ - ip{r). (6.4) 
Substituting Eq. into Eq. we obtain 
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dr(V'^(r)/(r)VXr)> = / drlV'tW + (pt(r)]/(r)[^(r) + ^(r)]. 



The right-hand side shows that there exists interference, as usual, it is described by the cross terms. 



V't(r)/(r)^(r)+^t(r)/(r)^(r). 

I 



(6.5) 



(6.6) 



In terminology of quantum optics, it represents the 
single-particle interference. This single-particle interfer- 
ence will vanish if the system goes into the normal phase 
where </?(r) = 0. A special case is the interference ap- 
pearing in the particle density where /(r) = 1, 



i/:^ {r)ip{r) + <y9t(r)i/)(r), 



(6.7) 



this form of interference is very familiar in optics, Eq. 
(|6.6() being the general case. 

In particular, an abnormal form of single-particle in- 
terference can appear in the electron-pair amplitude. 



(^X (r)^T (r)) = [V'i (r) + n (r)] [^T (r) + (r)] , (6.8) 
it is described by the abnormal cross terms 

V'i(r)</'TW + 'Pi(r)V'T(r) + 'Pi(r)</'TW> (6-9) 



which have no counterparts in optics, in contrast to Eq. 
(|6.7|) . Observe Eq. (|3.3I) . one recognizes that the field 
93 (r) contains two parts, 



where (p^^\r) represents the creation part. 



(6.10) 



(^(+)(r) = 



4(r) 



(6.11) 



with 0j (r) and (/)| (r) being connected with the operators 

cjj.| and c^kj, respectively; and (p^~\r) the annihilation 
part, which is connected with the operators Ck| and c_ki- 
This separation of two parts leads us to 



V'i(r)¥>T(i-) + ¥'i(r)V'T(r) + Vi(i-)^T(i-)=^i(r)0l(r)-K(/.|(r)V'T(r) + (pi > {v)cj)\{r) + cj)]{r)ip\ >{r). (6.12) 

Clearly, there exists interference on the right-hand side, this interference induce the abnormal interference on the left- 
hand side. The abnormal interference has an important consequence, that is, the electrons arc formed into Cooper 
pairs in the superconducting phase. Evidently, the abnormal interference vanishes in the normal phase. 

Also, there exists the two-particle interference in the superconducting phase, which arises from the statistical average 
of a two-body operator, 

^jj drdr'(V'^(r)V'^(r')w(r - r')V'(r')V'(r)) = ^ // drdr't;(r - r')i^^r)ii^r')i^{r')i^{r), (6.13) 



where v{r — r') stands for the two-body operator in 
Schrodinger picture. 

In the sense of Eqs. H6.5|l and (|6.13|) . the supercon- 
ducting phase can be identified as a form of interference 
pattern. Physically, the onset of an interference pattern 
will make a system more structured, and thus increases 
its degree of order and decreases its entropy. That is why 
the superconducting phase has a higher degree of order 
and a less amount of entropy than the normal phase. 

From the above discussions, it follows that there will 
appear interference in the ordered phase whereas there 
will not in the normal (disordered) phase. Obviously, the 
same conclusion holds for other phase transitions. There- 



fore, a phase transition amounts to the onset of a form 
of interference. As well-known, interference is a charac- 
teristic property of wave, we thus conclude that phase 
transitions originate physically from the wave nature of 
matter. 

Physically, any system will go in its disordered phase 
at sufficiently high temperatures, or as T — > -l-oo. Be- 
cause order parameter equals zero in disordered phase, 
the extended ensemble theory should ensure that there 
exists a zero solution of order parameter, at least at high 
temperatures. So it does, indeed. 

Upon a transformation, Eq. H2.6|) can be expressed as 
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S{(b, (3) = -Tr(ln(p(e-^^(^'^)ii' (c) e'"^*'^^)) p{H{c))^ . (6.14) 

Considering the expansion, 

^-^Di^,c)JJ ^^D(4..,c) = H (c) - i[D{^, c), H (c)] + c) , c), H (c)]] + • • • , (6.15) 

we obtain 

S{4>,[3) = -Tr(\n{p{H {c)))p{H{c))) + pt^T,{[D{4>,c),[D{4>,c),H {c)MH{c))) + ■ ■ ■ , (6.16) 

the hnear term of cf) vanishes. 

Apart from an irrelevant term, this equation shows that the powers of (f) are equal to or higher than the 2nd in the 
expansion of S{(f>,f3). As a result of 6S = 0, the powers of (j) are equal to or higher than the 1st in the equation of 
order parameter. It implies that there always exists a zero solution of (j) at any temperature, i.e., the trivial solution. 
In fact, this solution is a stable solution as T ^ +oo, 

AS ^ S{dcl>, (3) - S{0, (3) ^ -Tr(^ln{p{H{c))) [p{H'{Scf,,c)) - p{H{c))]) =0, T ^ +cx3. (6.17) 

I 



which satisfies the criterion of Eq. I|2.5|l . Of course, this 
conclusion is also valid for other phase transitions. 

To summarize, there exists a zero solution of order pa- 
rameter at any temperature in the extended ensemble 
theory, this solution becomes stable as T +oo. There- 
fore, a system always stays in its normal phase at suffi- 
ciently high temperatures, this high-temperature phase is 
disordered and structureless. As temperature decreases, 
there can arise spontaneously the quantum interference 
of matter waves, which makes the system structured and 
transform into its ordered phase. That is the physical 
picture for phase transitions within the framework of the 
extended ensemble theory. 

Remark: In Ref. [HO], Yang suggested an important 
concept, i.e., off-diagonal long-range order (ODLRO), 
to describe quantum phase transitions. Physically, 
ODLRO is a generalization of the concept of Cop- 
per pairing. With this generalized concept, he devel- 
oped many beautiful theorems, which are character- 
istics of quantum phases. Regrettably, in that pa- 
per, there is no discussion of the properties of the 
Hamiltonian that is needed to ensure the existence of 
ODLRO at low temperatures. Evidently, there exists 
no ODLRO if a system Hamiltonian is in its symmet- 
ric representation. The requirement can be satisfied if 
the system Hamiltonian takes its asymmetric, or broken- 



symmetry, representation. Therefore, those theorems in 
Ref. 50] still hold in the extended ensemble theory. 
As a matter of fact, the ODLRO 's of the reduced den- 
sity matrices pi and p2 Mh '3-g-> Tr(cicrp(ff'((A, c))cj^) 
and TT{ckcrCiap{H'{(j),c))Cj^cl^), result physically from 
the single-particle and two-particle interferences, respec- 
tively. That is to say, off-diagonal long-range order is 
a physical manifestation of the quantum interference of 
matter waves. 



7. INTERACTING BOSE SYSTEMS 

In this section, we focus our attention on the interact- 
ing Bose systems. We shall first study the weakly in- 
teracting Bose gas, and then the double-well potential 
systems where structural phase transitions, Goldstone 
bosons, and Higgs mechanism are concerned. 



7.1. Weakly Interacting Bose Gas 

As usual, the Hamiltonian for an interacting Bose gas 
reads. 



dr 



2m 



drdrVHr)V'Hr')t^(r " r')V'(i-')V'(r), 



(7.1) 



where ^/'(r) represents the boson-field operator, and v{r — r') the interaction. Expressed in terms of ip{r) and tp^{r), 
the phase-transition operator for BEC becomes. 



D[r],ip] = / dr [77^(r)V'(r) -H77(r)'(/'^(r)] 



(7.2) 
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With H[4>] and D [rj, tp], the entropy of the system can be obtained, 
5[77,/?] =-Tr(ln(p(i/ [7/^])) 



/? / dr 



— V77^(r) • V?7(r) - ^77^(r)77(r) 



/ / drdr'?7^(r)?7(r)t;(r — r')'0t(r')-0(r') 



+ j j drdv''q\v)ri{v')v{v - v')Tp^ {v')iI){y) + ^ JJ drdr'r]\r)ri{r)v{r - r')ri'^ {r')r]{r')^ , 



(7.3) 



where F denotes the statistical average of the operator F with respect to H[ip]. Its variation yields the equation of 
order parameter, 



2m 



/i ) 77(r)+ / c?r'7/(r)u(r— r')7/;t(r')'0(r')+ / dr'7/(r')u(r — r')'i/;t(r')'!/)(r)+ / dr'77(r)i;(r— r')77' (r')r;(r') = 0. 

(7.4) 

That is a generalized Ginzburg-Landau or Gross-Pitaevskii equation. In the simple case where u(r — r') — g6{r — r'), 
it reduces to 



2m 



VVr)+ 25Vnr)V'(r) - /X 77(r) + g ^(r)^ r7(r) = 0, 



(7.5) 



r 



which is the standard Ginzburg-Landau [5l| or Gross- interaction between two particles must be repulsive for an 

Pitaevskii [s^ equation. In addition to Eq. H7.4|) . interacting Bose system to be stable, i.e., v{r — r') > 0. 

one also needs Otherwise, v{0) < 0, the system will collapse. 

/Suppose that the repulsive interaction is weak, we can 

dr77^(r)?7(r) -I- drip'^ {r)^p{r) — N, (7.6) then treat it as a perturbation. As a technique of pertur- 

J bation. we shall emnlov GF. and adont the aelf-consisteni 



which is the equation of chemical potential. Eq. (|7.4|) 
and Eq. (|7.6|l constitute the two basic equations for an 
interacting Bose system. 

For simplicity, let us first consider the homogeneous 
solution. 



bation, we shall employ GF, and adopt the self-consistent 
Hartree-Fock approximation as in Ref. js^l , 



t/(k,ia;„) = 
KE* (k) = 



1 



iun - h-^ [e (k) - ^ - fiS* (k)] 
Hk' 

^[^-(0)+^.(k-k')] 



77(r) = 



(7.7) 



where ^ is independent of r. Accordingly, Eqs. (|7.3|l and 
(|7.4|) get simplified as 



{2n 
1 



^e*-""g(k',*c.„) 



(7.12) 



(7.13) 



S [v, (3] ^-Tr {In {p{H m)p{H W)) 
+ (3 j drh/ie^^+i«(0)(^t^)2] 

-M + v{Q)\^\\ = 0. 



where fi-S* (k) represents the proper self-energy. By use 
of Eq. 1)7. 10(1 . t/(k, iijjn) can be expressed as 



(7.8) 
(7.9) 



t/(k,ia;„) 



-fi-M£(k)-M]' 



(7.14) 



where 



where 



^ = /i 



dk 

(2^ 



b(0) + «(~k)] 



t dk' 

g'(k) = e (k) + / — 3 [t;(k - k') - ^;(-k')] 
J \2it) 



--j5:e--^(k,..„). 



drw(r)e 



-ik-r 



(7.10) 
(7.11) 



i-)^e-""t;(k',*c.„). (7.15) 



f3h 



Eq. H7.14|l shows that e (k) and ^ are the renormalized 
energy and chemical potential, respectively. 

With the help of Eqs. ^TT^ and (TTT^ . Eq. can 
Q{k,iLjJn) being the single-particle Green's function de- be reduced as 
fined with respect to H[ip]. 

Eq. (|7.8|l has the typical form of landau free energy 
[2l| . From Eqs. (ESJ, (ESJl and ^T^, it follows that the 



)=n, (7.16) 
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where n = N/V is the particle density of the system. 

When /I < 0, Eq. H7.9|l has the only solution ^ = 0. 
As shown by Eq. (|7.8(l . it is a stable solution. In this 
case, the system is in its normal phase. The renormalized 
chemical potential /I is given by Eq. H7.16|l . 



+ 00 



(7.17) 



where V denotes the principal value introduced in Ap- 
pendix^ and N{uj, Jl) the density of states, 



dk 

(2^ 



(7.18) 



If ^ > 0, Eqs. H7.8I) and (|7.9|) show that the system will 
transform into a condensation phase where 



v{oy 



Substitution of it into Eq. (|7.1t)|) yields 



+00 2 



= n, 



(7.19) 



(7.20) 



which gives the renormalized chemical potential /i at T < 
Tc, where Tc is the critical temperature for BEC, 



+ 00 



dtj7V(w,0)- 



1 



(7.21) 



that is to say, ^(T^) = 0. In sum, Eqs. l(7T7jl and (|7!^ 
yield together the renormalized chemical potential of the 
system at T > Tc and T < Tc, respectively. 

In order to deduce Jl and ^ in more detail, let us sup- 
pose further that v{r — r') is of short range. Following 
Ref. 0, we can approximate w(k) as 



v(k) = v{0) 



1- -(fca)2 



where 



J drr^v{r) 
Jdrv{r) 



(7.22) 



(7.23) 



Inserting Eq. iTT^ into Eq. (TTT^ and Eq. (TTTUI) 

leads us to 



2m* 



X V 



2v{0) 

+ 00 



v{0) 



nma m 



dLL>J^{LU, /i) 



3h? m 
(uj + Jl) 



where m* is the renormalized mass, 



1 

m* 



1 

m 



1 



v{0)nma'^ 
3fi2 



1 



(7.24) 



(7.25) 



(7.26) 



and Afiuj, the density of states defined by Eq. H7.18|l . 

AA(c.,u) = / 4^(2m*)3/2(.; + M)l/^ u;>~^ 
' [ 0, uj < —fi. 

(7.27) 

From Eqs. ^TT^, ifT?^ and (|T77|) . it follows that 



Tc v{0)nma'^ 



(7.28) 



where Tq stands for the critical temperature of the BEC 
happening in the corresponding ideal Bose gas, 



mkB VC(3/2) 



2/3 



(7.29) 



Eq. H7.28|l manifests that a repulsive interaction will 
lower the transition temperature of BEC, 



To To 3ft^ 



Also, we have by Eq. (frT7| . 



X 



1/2 



O /■+00 
1 ^3/2 / J 

0FC(3/2) io e-Ai/t-i 



(7.31) 



where t — T /Tc and ji = ji/ (ksTc). This equation gives 
the JiatT > Tc. I f T < Tc, Eq. I^T^ must be used in 
place of Eq. lf77r7|l . 



1 



X r 



V^C(3/2) VI + 7 
da;- 



ry 



3/2 



■ (x-uy)/t _ I ' 



(7.32) 



where y — fi/ (u(O)n), and — v{0)n/ (kBTc). This equa- 
tion wiU yield both the Jl and |,^|^ = /i/i;(0) at T < T^. 

The numerical results of Jl and |^|^ are shown in Figs. 
4 and 5, respectively. We observe that the renormalized 
chemical potential Jl is a monotonically decreasing func- 
tion of temperature, at T = Tc itself continuous but its 
derivative not. This agrees with the T > Tc result given 
by Ref. [H^l and the references therein. When T < Tc, 
the chemical potential Jl is set zero in Ref. |54:], as Ein- 
stein 45,, did in the ideal Bose gas. In the extended 
ensemble theory, there is no mathematical limit on /i, it 
will change with temperature even if T < Tc so as to 
ensure the conservation of particles, as pointed out in 
Appendix 1X1 The behavior of |^| shown in Fig. 5 is just 
expected for an order parameter. 
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The internal energy of the system, E, can be derived 
conveniently from E = (H) + fiN, 

2 



e = 



V^C(3/2) 



+00 3/2 



(7.33) 



e = v 



2v V^C(3/2) Vl+7-7y 
a; 



1+7 



5/2 



X P / da; 
Jo e 



where 



, -w, — , r < Tc, (7.34) 



1 E 



Combination of Eqs. H7.31|l and (|7.33|) gives Cy, the 
specific heat at constant volume, 

Cv 1535/2(2) 9^3/2(2;) 



where 



4 53/2(2) 451/2(2) 



T > Tc, (7.36) 



1 /■+°° x"^""^ 
dx- 



r(f^) Jo 2 ^e^ — 1 



(7.37) 
(7.38) 



This result is the same in form as that for the ideal Bose 
gas, which is familiar in standard books on quantum sta- 
tistical mechanics. If T < Tc, Cv must be calculated 
from Eqs. lf73^ and (173^ . 



The numerical results of Cy are shown in Fig. 6. The 
shape of Cy shows that the BEC occurring in the inter- 
acting Bose gas is a "A" -transition, as is naturally ex- 
pected. In contrast to the ideal Bose gas, the weakly 
interacting Bose gas manifests itself with two new fea- 
tures: First, Cv is discontinuous at T = Tc; and second, 
Cy is linear in temperature when T Tc- Evidently, 
this theoretical result for Cv meets the requirement of 
the third law of thermodynamics: Cy ^ as T ^ 0. 
The discontinuity aX T — is obviously due to the dis- 
continuity of the derivatives of Ji and |^| with respect 
to T. As to the linearity at low temperatures, it can be 
explained as follows. 



Observing that 



25 




FIG. 5: The order parameter for BEG versus temperature, where t — T/Tc, v = v{fS)nl (ksTc), and 7 = AT /To 
—v{Q)nmo? / (Sft'^). For the sake of contrast, the order parameter for the ideal Bose gas (i^ = 0, 7 = 0) is also plotted. 



r+ca ^1/2 
V / dx-, Tjr 

we obtain 



_ 2 



(7.40) 



r 



to the second power of t as t ^ 1. Inserting it into Eq. 
(1722), we have 



1 + 7 



V^C(3/2) Vl+7-7y 



3/2 



(7.41) 



it is an iterative equation oi y. At absolute zero (t = 0), 
it reduces to 



2/0 = 1- 



1 + 7 



3V^C(3/2) Vl+7-7yo 



3/2 



(^.yo)'/' . (7.42) 



This gives the zeroth order solution of y, which we de- 
noted by ?/o- Iterating of Eq. (|7.41|) to the second power 
of t leads to 



y^yo 



2C(2) 



1 + 7 



3/2 



V¥C(3/2) Vl + 7-7yo. 

(7.43) 

which shows that y, i.e., Jl and |^|^, decreases as in 
the low temperature region (T <C T^) , as can also be seen 
directly from Figs. 4 and 5. With the same procedure, 
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FIG. 6: The specific heat of a weakly interacting Bose gas versus temperature, where t = T/Tc, v — v{0)n/ (kBTc), and 
7 — at/To — —v{Q)nmo? / (3?i^) . For the sake of contrast, the specific heat of the ideal Bose gas [v = Q, ^ — Q) is also plotted. 



one finds 

f +00 



.3/2 



V I dx-, TTT 



(7.44) 

From Eqs. (|T^ . (|T^ . and l(T^ . it follows that e 
increases as T^. Therefore, the specific heat Cv will be 
linear in temperature T when T <^ Tc- 

If one expands Eqs. (|7.4l)|l and (|7.44l) further to the 
fourth power of t, he arrives at 



Cv = a(V^)T + b{V)T^, 



(7.45) 



where the expansion coefficients a{V) and 6(T^) are func- 
tions of V. Appendix IbI shows that this expansion does 
not depend on the approximations used here, it is valid 
provided that the interaction is weak. Therefore, it is a 
fundamental property of the weakly interacting Bose gas 
that the specific heat Cv vanishes linearly as T — 0. 

In connection with the A-transition occurring in liquid 
''He, obviously, the shape of Cv agrees well with that 
of the specific heat observed experimentally [s^. How- 
ever, the latter is Cp, the specific heat at constant pres- 
sure rather than at constant volume. As is well known 
from thermodynamics, they differ from each other by a 



temperature-dependent term. 



Cp=Cv+T 



dp 
df 



dV 

V ydTyp 



For solids at low temperatures. 



because 



Cp ~ Cv 



dV_ 
df 



(7.46) 



(7.47) 



0. 



That is why the Debye law for the phonons, the T 
law for the electron gases in normal metals, and the ex- 
ponential law for BCS superconductors are observed ex- 
perimentally in Cp albeit they are the laws of Cv at 
low temperatures. In other words, the theoretical result 
of Cv for a solid can be verified straightforwardly by 
the experiments. However, for a liquid, the situation is 
different, both the volume and pressure are much more 
sensitive to temperature. Therefore, the second term on 
the right-hand side of Eq. (|7.46|l can not be neglected 
any longer, the difference between Cp and Cv becomes 
important and should be taken into account. As a result. 
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to compare the theory with the experiment for a hquid, 
the experimental data for Cy or the theoretical result for 
Cp are needed. Regrettably, Cy is difficult to measure 
experimentally, and Cp is difficult to calculate theoret- 
ically. It is thus hard to compare the theoretical result 
with the experimental data. Fortunately, Cv tends to 
zero linearly as T ^ for a weakly interacting Bose gas, 
we can compare the theoretical result with the experi- 
mental data in the limit T 0. That can be seen clearly 
from the following analyses. 

According to the third law of thermodynamics, we have 



Sth = Sth{T, V) 



Cv{T\V) — , (7, 



49) 



where, as pointed in Sec. 13 Sth represents the ther- 
modynamical entropy. Substituting Eq. I|7.45|l into this 
equation, we obtain 



Sth ^ a{V)T + h{V)T\ 



(7.50) 



at low temperatures. Making use of the Maxwell rela- 
tions. 



\dTyv 
dV 
dT 



th 



dV 

dSth 
dP 



one has 
dP\ 

dV_ 
df 



a'{V)T+h'{V)T^, 



1, 



a'{V)T + ^b'{V)T^ 



p ydPJr 
Substituting them into Eq. (|7.4()|) leads us to 

Cp = a{V)T + b{V)T^, 



(7.51) 
(7.52) 

(7.53) 
(7.54) 

(7.55) 



where the terms with power higher than cube are omit- 
ted, and 



biV) = biV) [a'iV)f 



(7.56) 



Eq. (|7.55|) presents the form of Cp at low temperatures 
within the present theory. It contains a linear and a cubic 
term of T, the former includes only the contribution from 
Cv, and the latter includes both the contributions from 
Cv and the second term on the right-hand side of Eq. 
Heeding that Sth = at T = 0, we have 



mN 



dP 
dV 



(7.57) 



T=0 



where u denotes the velocity of sound. This implies that 
at low temperatures the velocity of sound is approximate 
to 



mN 



dP 
dV 



(7.58) 



That is to say, the factor {dV/dP)T on the right-hand 
side of Eq. H7.56|l represents physically the effect of 
sound. Therefore, the cubic term of Cp includes the 
contribution from the sound, which is in accordance with 
the viewpoint of Landau on quantum Bose liquid He II 
[HtI IH^ . This law has already been observed exper- 
imentally in liquid "^He at low temperatures [55l |. but it 
is not an intrinsic property of the BEC because sound 
exists in all fluids. The intrinsic property of the BEC is 
the linear term of Cp, it comes purely from Cy, and is 
characteristic of a weakly interacting Bose gas at T ^ Tc, 
as mentioned above. There is no such linear term in the 
Landau theory of quantum Bose liquid [s^ Is^ . To our 
knowledge, this linear behavior of Cp has not yet been 
reported experimentally for liquid ^He, it is thus a pre- 
diction of the present theory. Eq. H7.55|l indicates that 
the temperature T has to go much lower for Cp to show 
the linear behavior than to show the law. We believe 
that it would be observed if the measuring temperature 
is lowered enough, and if liquid '^He could be regarded as 
a weakly interacting Bose gas and the A-transition were 
a Bose-Einstein condensation. 

Remark: In general, |^|^ can not be interpreted as the 
density of condensed particles. This can be seen from Eq. 
(ESS, it shows that ?/ > 1 at T = OK, that is. 



(7.59) 



at absolute zero. Only in the case of the ideal Bose gas, 
1^1^ can not be greater than the density of particles, i.e., 



(7.60) 



at any temperature. Physically, f is the order parame- 
ter and internal spontaneous field of the system, there 
is no a priori reason why it must be interpreted using 
the the particle density. In fact, any statistical average 
of observable, including the particle density, will depend 
explicitly on ^ when T < Tc- 

Now, let us discuss the correlation of the BEC to the 
superfluidity. As usual, we set 



(7.61) 



where p{r) = \r]{r)\ and (p{r) = arg(?7(r)). With this, 
we have 



j(r) 



ih 
2m 



((V-V>t(r)v,(r'))|r'. 



— p(r)V(^(r) 



(7.62) 



obviously, j(r) is the supercurrent density. If (^(r), the 
phase of the BEC order parameter, is independent of r, 
there is no supercurrent, j(r) = 0; otherwise, j(r) ^ 0. 
The supercurrent velocity v is given by 



v(r) — — V(y3(r). 

m 



(7.63) 
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Since the order parameter ?7(r) is a single- valued function 
of the position r, we have 



/v(r) ■dr = l — , le Z. 
J m 



(7.64) 



That is to say, there can exist quantized vortices in the 
condensed phase if (p{r) depends on r. We remark that 
there exists a critical velocity Vc, the velocity of the su- 
percurrent can not be greater than it, i.e., v < Vc- That 
is because the increase of the supercurrent velocity is un- 
favorable to the decrease of the entropy of the system; 
the larger the v is, the greater the entropy will be, which 
can be easily seen from the first integral over r on the 
right-hand side of Eq. (|7.3|l with, e.g., ri{r) — ^e^^'^ 
where p = constant. 

Finally, we shall prove that there can not exist any 
supercurrent or quantized vortex in the ideal Bose gas. 

For the ideal Bose gas, v{v — r') = 0, Eq. I|7.4|l reduces 

to 



2m 



(7.65) 



obviously, it is the eigenvalue equation of the kinetic en- 
ergy operator. This equation can be rewritten as 

/ dr\—ihVri{r)P ■ [—ihVr]{r)] — IJ. [ dr ?7^(r)?7(r). 

2m J J 

(7.66) 

It indicates that, for any nontrivial eigenfunction, the 
chemical potential ^ can not be negative, i.e., > 0. 
Therefore, there can only exist the trivial solution, ry(r) = 
0, when /i < 0. The nontrivial solution, r]{r) ^ 0, can 
exist only when /i > for the ideal Bose gas. 

To discuss the stabilities of those solutions of Eq. 
(|7.65l) . one needs to investigate whether AS > 0. From 
Eq. (lESJ, we find 



AS ^ 13 dr 



2m 



VStji^r) ■ V(5?7(r) - n6r)^r)Sr]{r) 



(7.67) 

where Sri{r) represents the variation of the order param- 
eter 7y(r) from any one of the solutions of Eq. H7.65|l . 
Evidently, AS > when p < 0, this means that the sys- 
tem is stable when fj, < 0. Since 77(r) = when < and 
r](r) = constant when /i = 0, the normal phase and the 
homogeneous condensed phase are both stable, which is 
in accordance with the results of Sec. 15.1.21 As to the 
case of ^ > 0, the solution of Eq. (|7.65|) is inhomoge- 
neous, which is just the case we are now interested in 
because an inhomogeneous solution can produce super- 
current and quantized vortices. To determine whether 
there can exist any supercurrent and quantized vortex 
in the condensed phase, we first reformulate the above 
equation as follows. 



AS ^ f3 dr 



' ^ ' ^ 2m 



(5ry(r) — fiSi]' {r)Sri{r 



(7.68) 



and then examine the eigenvalue problem, 

V^(5?7(r) = ESr](r), 

2m 

(5?7(r)|, 



r^+oo 



0. 



(7.69) 
(7.70) 

Eq. H7.70|l is the boundary condition obeyed by the vari- 
ation Sri{r). Obviously, this eigenvalue problem always 
has solutions when E > 0, e.g.. 



2m 



(7.71) 



where ji{z) is the spherical Bessel function. Accompany- 
ing those solutions, AS can be reexpressed as 



clearly, 



AS^P driE - n)Sr]^r)Sr]ir) 



> 0, E> p 
AS = 0, E = pL 
AS <0, E <p. 



(7.72) 



(7.73) 



It shows that any solution of Eq. (|7.65(l for /i > be- 
longs to the saddle points of the system entropy. In other 
words, any inhomogeneous condensed phase is physically 
instable, needless to say, the supercurrent and quantized 
vortex. This proves that there can not exist any super- 
current or quantized vortex in the ideal Bose gas. 

By the way, the above analyses demonstrate that the 
condensed phase of the ideal Bose gas must be homoge- 
neous. There is no inhomogeneous condensation in the 
ideal Bose gas. This result is very natural, it confirms 
the Einstein's deep physical insight on BEG, again. 



T.2. Double-well Potential and the BEC in 
Configuration Space 

From Sec. 15.31 it is learned that an ideal phonon gas 
cannot produce BEC. However, the ideal phonon gas 
is simply a harmonic approximation to an actual solid. 
Generally, the interatomic interaction is anharmonic. It 
thus raises the question as to whether an anharmonic sys- 
tem can produce BEC. We intend to discuss the question 
in this subsection. 

To that end, it is helpful to make a representation 
transformation to the formulation used in Sec. 15.31 
There, what we used is the Fock space where all the 
observables are expressed in terms of creation and an- 
nihilation operators. In the following, we prefer to use 
the phase space where all the observables are expressed 
in terms of generalized coordinates and momenta. As 
is well known, the two spaces are equivalent and can be 
transformed into each other according to such a rule |59j| , 



qi = 



y 2miUJi ^ 



(7.74) 
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and 



Commonly, a system Hamiltonian has the form, 



y/2/niihuJ. 



-Pi, 



(7.75) 



where Qi and Pi are the zth pair of coordinate and momen- 
tum; a.i and aj the corresponding annihilation and cre- 
ation operators; rrii the effective mass; and Ui > the zth 
parameter, which can generally take any positive value, 
and particularly the natural frequency if qi and pi consti- 
tute a harmonic oscillator. Along with Eqs. H7.74|l and 
(|7.75l) , the main quantities concerned are transformed as 
follows, 



F{a) 
H{a) 
D(e,a) 



Fiq,p), 
Hiq,p), 



where 



i 

D{r],C,Q,P) = + QPi): 



(el- 6). 



h_ 



(7.76) 
(7.77) 
(7.78) 
(7.79) 



(7.80) 
(7.81) 

(7.82) 



Here, ^ represents the order parameter for BEC in Fock 
space, rii and Q the corresponding order parameters in 
phase space. 

Under the action of D{ri, q,p), the statistical aver- 
ages of qi and Pi will change with temperature as 



(<Z,;> = -C^+Tv{qip{H{q,p))) 
{Pi) =Vr+ T'r{p,p{H{q,p))). 



(7.83) 



If the system produces a BEC, i.e., rji and/or Q change 
from zero to nonzero with the decreasing of temperature, 
the average positions of the particles will redistribute in 
phase space. In this sense, the BEC is said to be the 
phase-space BEC, it changes the distribution of the sys- 
tem in phase space. If, however, the BEC only causes rji 
(or Q) to change from zero to nonzero, the particles will 
redistribute merely in momentum space (or configura- 
tion space), the subspace of phase space. In such a case, 
the BEC is said to be the momentum-space BEC (or the 
configuration-space BEC), it changes only the distribu- 
tion of the system in momentum space (or configuration 
space). Obviously, once a BEC changes the distribution 
of a system in configuration space (Ci ^ 0, 77^ = or 
not), it induces a structural phase transition (SPT) '6O]. 
In this connection, it can be said that a SPT is just an 
instance of BEC. 



'^u{q,) + -'^v{q^ 



-q,), (7.84) 



where u {qi) and v{qi — qj) stand for the single- and two- 
particle potentials, respectively. One can easily verify, 
using Eqs. ifT^ . ((775|) and ^TJB: that r]i must 

be zero at any temperature, i.e., there cannot arise the 
momentum-space BEC. In consequence, the system with 
such a form of Hamiltonian can, at most, produce a 
configuration-space BEC. Therefore, a SPT is commonly 
a configuration-space BEC. In the following, we shall con- 
fine our attention within this kind of SPT. The phase- 
transition operator D(r], C, q,p) can now be simplified as 



(7.85) 



for rji vanishes forever. 

In the field of SPT, the Hamiltonian of Eq. (|7:M|) is 
made concrete as l6(l. 



i i 



1 2 2 , 1 4 
2^ 9» + -^9<lt 



[7.. 



where /i^, g, and A are constants, and (ij) denotes the 
nearest neighbor sites. The system gets anharmonic if 
fi^ < 01 g ^ 0. The terms included by the second sum 
represent the single-particle potential that arises from an 
underlying sublattice of atoms, which do not participate 
actively in the phase transition, and the one included 
by the third sum represents the elastic coupling between 
nearest neighbors, which is described by a harmonic po- 
tential with the coupling constant A > 0. Obviously, 
H{q,p) has the discrete symmetry. 



it leads us to 
where 



s ■ qie ' = -qi, 

(7.87) 



9l = 0, 



Fiq,p)^Tr{F{q,p)piHiq,p)y 



[7.. 



(7.89) 



This discrete symmetry will break down if a SPT takes 
place. 

Let us begin with the simple case of A = 0. Making 
use of Eq. (|7.85|) , we have 
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3.9 (q 



If 



(7.90) 



r 



It shows that the entropy S{C,, (3) has the typical form of 
Landau free energy |2l| . From Eqs. 
(|2.5|) . one can draw three conclusions. 



GIOI, EH), and 



good basic part and is thus unfavorable for doing ap- 
proximations. Complying with quantum field theory, the 
transformation can be achieved as follows, 



1. The case of 5 < 0. The system is instable and will 
collapse because the entropy is unbounded from be- 
low with respect to the variation of Q^. 

2. The case of 5 — 0. (a) If p^ > 0, the trivial solution 

= is the only minimizer of S{(,/3). That is to 
say, a system consisting of harmonic oscillators, or 
rather an ideal phonon gas, can not produce SPT, it 
will stay in the normal phase forever, which is in ac- 
cordance with the result of Sec. 15.31 (b) If = 0, 

can take a temperature-independent value, it 
represents, as indicated by Eqs. H7.83|l and H7.88|l . 
the equilibrium position of qi, i.e., (qi) = —Q. As 
the initial equilibrium position of qi at f3 = 0, —Q 
can be chosen as zero. Therefore, a system com- 
posed of free distinguishable particles can not pro- 
duce SPT, as is just expected, (c) If < 0, Ci = 
is the only solution, but it is a maximizer of 5(C, /?), 
the system is instable and will collapse. This result 
is quite natural since all particles go in a reversed 
harmonic potential. 

3. The case of g > 0. (a) If > 0, G = is the only 
solution, and it is a minimizer of S'(C,/9). There 
is no phase transition in the system, the potential 
is a single well, (b) If /i^ < 0, the potential gets 
double welled, and it must be determined by the 
position fluctuation q'^—'qf' whether the system can 
produce a phase transition. Evidently, a transition 
can occur if the fluctuation is such a decreasing 
function of temperature that /i^ -I- ig{q^ — qf) can 
change in sign, from positive to negative. If the 

fluctuation is so strong that ^J-"^ + 'ig[{qf) — ijii) ] > 
up to zero temperature, the system has to stay in 
the normal phase forever. 

From the above discussions, it follows that a double- 
well potential where < and g > is the only possible 
case for the system to produce a structural phase transi- 
tion under the condition of A = 0. One can easily verify 
that this conclusion is still valid for A > 0. 

We are now confronted with performing the evalua- 
tion of the position fluctuation qf — ql^ so as to ascertain 
whether a phase transition can occur or not if < 
and g > Q. To this end, we need to perform a represen- 
tation transformation to the Hamiltonian of Eq. (|7.8t)|) 
because, as well-known, this form of Hamiltonian has no 




= e ^ qie ' = qi 

TTi = e ^ PiC ' ^Pi 



where 



(7.91) 



(7.92) 



In the sense of classical mechanics, the above transfor- 
mation amounts to performing a translation so that the 
new coordinate (j)i is referenced from qi = — i/, one of 
the two minimizers of the double-well potential. Along 
with Eq. (|7.91|) . the new representation H{(j),TT) of the 
Hamiltonian becomes 



iY.'^Pj -iE'^Pj 
— e ^ H(q,p)e ^ 

4 



1 4^ 



(7.93) 



The first sum on the right-hand side represents harmonic 
oscillators with natural frequencies uji — ^ jm^ it 
constitutes a basic Hamiltonian which is favorable for 
approximation handling, the second the interaction, and 
the third an unimportant constant. 
With Hid), tt), one finds 



~2 



(7.94) 



where F stands for the statistical average with respect to 

H(0,^), 



Tr(F(g,p)p(i/(g,p))), 



(7.95) 



where H{q,p) = ff((/>, tt). According to Eq. H7.94|l . we 

now just need to evaluate the fluctuation qf — qi^ , it can 
be worked out by the retarded Green's function {{qi \qi))u: 
which is defined with respect to H{q,p) [s^l- The Green's 
function satisfies the following equation of motion, 



f^{{qi\qi))Lo = + {{[qi,H{q,p)]\qj))^ 

= —{{pi\qi))u- 

m 



(7.96) 
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Also, 

fy^{(jpMi))u> {[Pi.qi]) + {{[Pi,H{q,p)]\q,))^ 

= -ih[[l - 2^l''{{qM^))u + 3gi^{{qfh)).j 
+ 9{{qfh))^]- (7.97) 

We would like to truncate the GF chain with the follow- 
ing decoupling approximation, 

{{qlh))u. = mq^Wi))^ = -i^((9.k.))c., (7.98) 

mq^))^ = qt{{q^\q^))^ = [(^ - ft') + ^']((ftlft))^, 

(7.99) 

where we have used the result — qi + u = which 
follows from Eq. H7.88|l . This approximation leads us to 



((ftlft))c^ = 



1 1 



where 



m 



1/2 



(7.100) 



(7.101) 



represents the renormalized frequency. According to the 
fluctuation-dissipation theorem Isslj , the fluctuation q^ — 
qf' can be obtained by 



q'-qf^V 



it gives rise to 



ql -q^^ = 



+ 00 



duj- 



1 



Im((ftlft))a;+»0 



1 



1 



(7.102) 



(7.103) 



Eq. (|7.1U3|) and Eq. (|7.1U1|I constitute a self-consistent 
approximation to the fluctuation q^ —qf'. Evidently, the 
first term on the right-hand side of Eq. H7.103|) represents 
the zero-point fluctuation. Eq. (|7.103|l shows that the 
fluctuation gets weaker and weaker when temperature 
goes lower and lower. At zero temperature, one has 



ql - q> 



(7.104) 



and at high temperatures, 



From Eq. H7.106|) and Eq. ifT^^ . it follows that a SPT 
will occur at r„, 



-tL. 

9 



3h 



1/2 



T>Tc 

T < Tr. 



(7.107) 

Specifically, Q = —p? j g at zero temperature if the 
zero-point fluctuation is neglected, i.e., the order param- 
eter will equal one of the two minimizers of the double- 
well potential at zero temperature when the zero-point 
fluctuation is left out of consideration. 

Combining Eqs. (|7.1()7|l . ((71^ . and ifTl^ . we have 



0, 



li_ _ _3h_ ( 1 

g inuJs I 2 



1/2 



T>T, 

T<n. 



2 ' JFfi^s _i ^ ^ 

(7.108) 

It shows that for each site the distortion (g^) can take 
independently either of the two values at T < Tc, there 
is no correlation between the distortions of different sites, 
therefore, the new phase below Tc is a structural glass. 
Physically, that is because we have discarded the coupling 
between nearest neighbors (A = 0). If it is taken into 
account (A > 0), one can easily verify that the new phase 
below the transition temperature will be ferrodistortive: 



(qi) = (qj) 



(7.109) 



where i and j denote any two different sites. This result 
is in agreement with the classical mean-field theory [s^l . 

In conclusion, we find that the double-well poten- 
tial system described by the Hamiltonian of Eq. (|7.86() 
can produce a structural phase transition. Mechani- 
cally, that is because the position fiuctuations of particles 
will decrease with temperature. In the high-temperature 
regime, the fluctuations are relatively strong, the system 
can only stay in its normal phase. The fluctuations will 
get weaker and weaker when temperature goes lower and 
lower, and finally the system undergoes a transition at 
a finite temperature and then turns into a new ordered 
phase in the low-temperature regime. We learn from here 
that the position fluctuations of particles play an impor- 
tant role in structural phase transitions. Also, we see 
that configuration space is very convenient for describing 
SPT, in contrast to the Fock space. 



7.3. Goldstone Bosons, Ginzburg-Landau 
Equations, and Higgs Mechanism 



q!~qf = 



1 



gp 



0{T 



l/2^ 



(7.105) 



They imply that there exists a finite temperature Tc such 
that 



A*' + 3g(gf-gl') <0, T < T,, 



(7.106) 



Obviously, the preceding theory for the discrete case 
can be extended straightforwardly to the continuum case. 
Let us consider the so-called 0(A^)-symmetric vector 
model (61. .62]. 



H [V', tt] = / dx 



+ ^^V + f(^W 



(7.110) 
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where < and g > are the two constant param- 
eters for a double- well potential, ij^i^) is usually called 
Higgs field, it represents a real- valued vector field with 
N components, and 7r(x) the corresponding momentum, 

[M^),T^ji^')]=^6,jS{^-^'). (7.11f) 

As is well known, this double- well Hamiltonian is the t^- 
ical model for SSB in quantum field theory 61. .62,. ,.63|. 
The arguments given in the preceding subsection are still 
valid. They reveal that the physical reason why the 



0{N) symmetry will break down is that the fluctuation 
of the field ip{:x.) decreases monotonically with tempera- 
ture, and that this SSB belongs essentially to the EEC 
happening in configuration space. 

In the condensation phase, the system Hamiltonian 
takes the form, 

= e'^I^'"li/[V',7r]e"^^['^'"l, (7.II2) 

where <y3(x) is the order parameter for Higgs field, and 
D [</?, tt] the phase-transition operator. 



i:»[V5,7r] =y dx (^^(x)7r(x). (7.113) 
In more detail, Eq. (|7.112|l can be written as 
H'[^,^,7r]= I dx |l^t^+ lv(^t+<^t).v(^ + ^) + !^(V-t+v5t)(V'-t-<^) + |[(^t + ^t)(^ + <^)]2l (7,114) 



At absolute zero, if one neglects the zero-point fluctu- is the minimizer of the double-well potential. Without 
ation of the field i-e., '4'^'>P ~ 0, then he has the stable loss of generality, we shaU take 
solution, 

'^^'^=^'' (^-115) ^ hO,...,0]. (7.117) 

where 



jy = yj-rri^jg (7.116) Substituting it into Eq. I|7.114|) leads us to 



dx 



11 v^' 
I 



(7.118) 



That is the zero-temperature representation of the sys- 
tem Hamiltonian, it shows that the first component field 
gets a mass of — 2m^, the others are massless. 
As well-known, those massless particles are the so-called 
Goldstone bosons From here, one can see how Gold- 
stone bosons are produced in a natural way by sponta- 
neous symmetry breaking within the framework of the 
extended ensemble theory. Of course, he should also rec- 
ognize that the zero-point fluctuation can impose a little 
influence on the Goldstone bosons. 



At last, it should be emphasized that, in the 0{N)- 
symmetric vector model, Goldstone bosons originate 
physically from the Bosc-Einstcin condensation of Higgs 
field. 

If a gauge field is coupled with the Higgs field tp{:xi), 
what will happen to the gauge field? can it produce BEG 
as the Higgs field iIj{x.)7 To clarify this problem, let us 
consider, for simplicity, the coupling of a complex-valued 
(or two-component) Higgs field ip{'x.) and an electromag- 
netic field A(x), 



H [■(/i, tt; A, E] 



dx 



(V -t- ieA) • (V ■ 



ieA) tp 

r 



\ (V X A)^ 



(7.119) 



where e is the charge of the Higgs field '0(x), and E(x) the electric field strength which plays the role of the canoni- 
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cal momentum corresponding to the canonical coordinate 
A(x), 

[A,(x), -E,(x')] = ^%(5(x - x'). (7.120) 

Here, the temporal gauge has been used to perform a 
canonical quantization to the electromagnetic field, an 
Abelian gauge field. The reason for choosing this gauge 
is that it can be easily applied to quantize non- Abelian 
gauge fields, and transformed into other gauges through 
the path-integral formalism introduced by Faddeev and 



Popov m. 

For this coupled system, BEC can be explored through 
the phase-transition operator, 

D [V3,7r; A,E] = J dx {tpn + tpK^ - A ■ F,) , (7.121) 

where ip{'x.) and A(x) are the order parameters for the 
Higgs field and gauge field A(x), respectively. With 
D [if, tt; a, E], one can obtain the entropy of the system. 



S [if, A,/3] = -Tr( ln(p(iJ [ifj, tt; A, E]))p(iJ [ifj, tt; A, E])) + / dx (V + ieA) (^^ • (V - ieA) ip 



+ (m^ + e^A • A -I- AgTp^il;)(p^ ip + g ((/sV) + e^^t^jA^ + 2 ^ ^) 



(7.122) 



which shows that the two order parameters <p and A are coupled together. Eq. H7.122|l reminds us of the Landau free 
energy for superconductivity |5lll56l|. indeed, they are both similar in form. 
Variation of S with respect to Lp and A yields 



- (V - ieA)^ (p + (m^ + e^A • A + '\:g^^ip)(p + 2gip^ tptp = 
V X V X A = -ie [ip^Vip - ipVp'^) - 2e^((^V + i^H)^- 

I 



0, 



(7.123) 
(7.124) 



This form of two coupled nonlinear equations is known 
as Ginzburg-Landau equations [5ll l56j | . Observing that 

(B) = -V X A, (7.125) 

where B denotes the magnetic field strength, one can 
specify 

V • A = 0, (7.126) 

which indicates that there are only two independent order 
parameters for the electromagnetic field. Eq. (|7.126|l can 
be regarded as the Coulomb gauge to the vector order 
parameter A. 

From Eqs. (|7.122trr.l2(j|) . one can easily verify that 
there is no BEC for the electromagnetic field A if it is 
not coupled with the Higgs field V', or it is coupled with 
ij) but the latter has not condensed yet {<p = 0). This 
also confirms the conclusion of Sec. 15.31 an ideal pho- 
ton gas, or a free electromagnetic filed, can not produce 
BEC. Once the Higgs field gets condensed ((^ 7^ 0), the 
electromagnetic field A must condense down (A ^ 0) 
simultaneously; otherwise, there can arise a linear term 
of 5A in 5S, the coupled system will go instable. This 
shows that the gauge field A will condense together with 
the Higgs fields when it is coupled with the latter. In 
such case, the transition is cooperative, one usually says 
that Lp is the primary order parameter, and A the sec- 
ondary one. 

With the help of Eq. ((71^ . Eq. (fTT^ can be 



reduced as 

V2A-2e^(V3V + 'i/^)A = -ie{p^V tp- ipV tp^). (7.127) 

Heeding that (A) = -A Eq. (|7.127|l is identical to the 
static Proca equation [63 except that the prefactor of A 
may be a function of r, the counterpart in Proca equa- 
tion just a constant. In general, the prefactor cannot be 
interpreted as the mass of the gauge field. However, once 
there is a constant occurring in the prefactor, a mass /i is 
generated to the gauge field, as can be seen more clearly 
in the following discussion. 

Owing to the nonlinearity, the Ginzburg-Landau equa- 
tions (|7.123|l and H7.124|) are rather complicated to han- 
dle. As usual, let us consider the zeroth-order approxi- 
mation of Eq. (TTT^ at T = OK mil, 

V3^(x)¥j(x) = (7.128) 

where the effect of A and the zero-point fiuctuations of 
the fields -0 and A are omitted. From Eqs. H7.127|) and 
(|7.128|l . it follows that 

V^{A) -e^v^{A) ^Q. (7.129) 

Since the prefactor e^v'^ is constant now, it can be inter- 
preted as the mass of the gauge field at zero temperature, 

;i2(r = 0) = 6^1/2. (7.130) 



34 



As well-know, this mass can account for the Meissner 
effect jl^l, i.e., the expulsion of a magnetic field from the 
interior of a superconducting material, with a London 
penetration length Al ~ "gS]. 

The above formalism is Ginzburg-Landau picture of 
mass producing, in this picture, the mass is represented 
by the equation of motion. There is another picture, 
which is due to Higgs ,69, j73,0|- In Higgs picture, the 
mass is represented directly by the system Hamiltonian 
itself. Since an equation of motion is in accordance with 
its Hamiltonian, both pictures are equivalent physically. 
To show Higgs picture, we aught to consider the zero- 
temperature representation Hq of the system Hamilto- 
nian, 



Ho = e'D^'fi'^'^'^^H tt; A, E] g-^^I'^^'^^^'El 
= ff[x,^;W,E], 



where 



(7.131) 

(7.132) 
(7.133) 



As is well known, within the temporal gauge, the elec- 
tromagnetic field still has gauge degrees of freedom, 
the Hamiltonian vr; A, E] remains invariant under 
any time-independent gauge transformation. As an- 
other representation of H[ifi,Tr; A,E], the Hamiltonian 
H [x, tt; W, E] is also gauge invariant, 

H [x, tt; W, E] = e'G[e,x.7r,E]H [x, tt; W, E] e-«G[e,x,^,E] ^ 

(7.134) 

where e M and 

G [9, X, TT, E] = J dx [~iee (ttx - x^t^^) + V6l • E]. 

(7.135) 

As Higgs [63, UM 'ii'l' '^^^ re- gauge the fields 
x(x) and W(x) by letting 0(x) equal to the phase of 
X(x), that is. 



x(x)=p(x)e'^W, p(x)e 



(7.136) 



so as to eliminate the Goldstone bosons represented by 
the phase of x(x). Upon such a choice of ^(x), the Hamil- 
tonian Hq can be expressed as 



iJo = / dx 



2 2 

m p 



gp 



(V + ieA)p ■ (V - ieA)p 



-E^ + -(V X A) 



(7.137) 



where 



(7.138) 
(7.139) 



A(x) = e'G[fl.x,-,E]^(^)g-^G[9,x.-,E]^ 

From Eqs. (|7.132|) and H7.121|) . it follows that 

X(x) =V^(x)-|-^(x). (7.140) 



Without loss of generality, we can, from Eq. H7.128|l . 
choose ip{x) as real. 



That leads to 



(^(x) = iy/V2. 



p(x) = (0(x) + i/) /V2, 



(7.141) 
(7.142) 



where 



(x) = V2e'^i'^^-'^^%i^)e-'^i''^'^^^l (7.143) 



Eq. (|7.142|l indicates that 0(x) is merely a real field. 
Substitution of Eq. I|7.142() into Eq. I|7.137|) gives rise 

to 



Ho 



dx 



^p' + ^(V + ieA)0-(V-ieA) 



m202 + 9^03 _^ + + i(V X Af 

^e^iy^A ■ A + e^vipA ■ A 

, 1 2 
— + 



where 

p(x) = [it(x) + i(x)]/V2, 
Po(x) =z[i^(x)-^(x)]/\/2. 

Paying attention to 

giG[e,x,'r.E] ^ X(x) + X^(x) ^ g-iG[e,x,7r,E] 



we have 



Kx),0(x')] = -^<5(x- 
bo(x),i/o] -0. 



(7.144) 

(7.145) 
(7.146) 



(X) V. 
(7.147) 

(7.148) 
(7.149) 



That is to say, the (/'(x) and p(x) become a new pair of 
canonical field operators whereas po(x) is just a constant. 
Physically, that is because po (x) corresponds to the elimi- 
nated degree of freedom (/'o(x) = —i (x(x) — X^(x)) /\/2- 

Eq. H7.144() forms Higgs picture of mass producing. 
Comparing it with Eq. (|7.129|l . one sees immediately 
that the mass of the gauge field in Higgs picture is iden- 
tical to that in Ginzburg-Landau picture, as is expected. 

The above theory for a complex field coupled with an 
Abelian gauge field can be easily generalized to the case of 
an iV-component field coupled with a non- Abelian gauge 
field, with the same conclusion that the gauge field will 
obtain a mass through spontaneous symmetry breaking. 
This manner of mass producing is the well-known Higgs 
mechanism, which was successively contributed by Lon- 
don l6a. Ginzburg and Landa u Isil . Anderson |72. IT^. 
HiggslilzSEIl, and Kibble [ziP 

Remark: Spontaneous symmetry breaking can occur 
without any fundamental Higgs field in superconductiv- 
ity, it is thus called " dynamical symmetry breaking" , 
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which imphes that Higgs fields are " normally" needed for 
symmetry breaking. From the viewpoint of the extended 
ensemble theory, that is a prejudice because supercon- 
ductivity, Goldstone bosons and Higgs mechanism share 
the same physical ground. This is in accordance with the 
standpoint of Huang expressed in his book [j^ . 



8. SUMMARY AND CONCLUSIONS 

So far, a possible extension of Gibbs ensemble theory 
has been postulated, as a personal review of the author, 
to enable a microscopic description to phase transitions 
and spontaneous symmetry breaking. The extension is 
founded on three hypotheses, which root, in physics, from 
Landau's ideas on phase transitions: order parameter, 
variational principle, representation transformation, and 
spontaneous symmetry breaking. In this sense, the ex- 
tended Gibbs ensemble theory can be viewed as a micro- 
scopic realization of the Landau phenomenological theory 
of phase transitions. 

Within the framework of the extended ensemble the- 
ory, a phase transition occurs according to the principle 
of least entropy, which manifests itself as an equation of 
motion of order parameter. This equation determines the 
evolution of order parameter with temperature, and thus 
controls the change in representation of system Hamilto- 
nian. Different phases correspond to different representa- 
tions, and vice versa. A system Hamiltonian will realize 
its symmetric representation in the disordered phase, and 
asymmetric one in the ordered phase. The change in sym- 
metry results from the change in representation. That is 
the Landau mechanism responsible for phase transitions 
and spontaneous symmetry breaking in the extended en- 
semble theory, it holds in the whole range of temperature, 
including the critical region. 

Physically, phase transitions originate from the wave 
nature of matter. A system always stays in its normal 
phase at sufficiently high temperatures, which is disor- 
dered and structureless. As temperature decreases, mat- 
ter waves will interfere automatically with one another. 
This interference makes the system structured and trans- 
form into its ordered phase. That is the physical picture 
for phase transitions in the extended ensemble theory. 

Also, the extended ensemble theory has been applied 
to typical quantum many-body systems, with the conclu- 
sions as follows. 

For the ideal Fermi gas, it can not produce supercon- 
ductivity, its normal phase is stable at any temperature. 

Negative-temperature and laser phases arise from the 
same mechanism as phase transitions, they are both in- 
stable, and will finally turn into the positive-temperature 
phase. That is the microscopic interpretation for nega- 
tive temperatures and laser phases. 

For the conventional weak-coupling low-Tc supercon- 
ductors, the mean-field solution due to Bardeen, Cooper 
and Schrieffer is derived anew, and proved to be stable 
within the extended ensemble theory. 



For the ideal Bose gas, it can produce Bose-Einstein 
condensation only in the thermodynamic limit, which 
agrees with Einstein's prediction. But that holds in the 
sense of Lebesgue integration rather than Riemann in- 
tegration. In general, the order parameter for BEG can 
not be interpreted as the number of condensed particles. 
Besides, there can not exist any supercurrent or vortex 
in the ideal Bose gas, the condensation is always homo- 
geneous. 

The ideal phonon gas can not produce Bose-Einstein 
condensation, neither the photon gas in a black body. 

It is not admissible to quantize Dirac field using Bose- 
Einstein statistics. Otherwise, the system is instable. 
That is a statistical rather than mechanical reason why 
Dirac field has to be quantized using Fermi-Dirac statis- 
tics. This conclusion is in accordance with the spin- 
statistics theorem in quantum field theory. 

A structural phase transition belongs physically to 
the Bose-Einstein condensation occurring in configura- 
tion space. For a double-well anharmonic system, it is 
instable at low temperatures, and will undergo a struc- 
tural phase transition at a finite temperature. Physically, 
that is because the position fluctuation is a monotonically 
decreasing function of temperature. 

For the 0(7V)-symmetric vector model, the 0{N) sym- 
metry will break down spontaneously, with the presence 
of Goldstone bosons. In essence, this SSB is a Bose- 
Einstein condensation. If the system is coupled with a 
gauge field, it can cause the latter to condense together 
with itself. The cooperative condensation can be de- 
scribed by the Ginzburg-Landau equations. After the 
condensation, the gauge field can obtain a mass, which 
is the so-called Higgs mechanism. 

For an interacting Bose gas, it is stable only if the in- 
teraction is repulsive. The BEG present in this system 
can be described by the generalized Ginzburg-Landau 
equation. If the interaction is weak, the BEG is a "A"- 
transition, and its transition temperature can be lowered 
by the repulsive interaction. As a characteristic prop- 
erty of this "A" -transition, the specific heat at constant 
volume Cv will vanish linearly as T —^ 0. 

If liquid ''He could be regarded as a weakly interacting 
Bose gas, and if the A-transition were a Bose-Einstein 
condensation, then its specific heat at constant pressure 
C p would show a law at low temperatures, which is in 
agreement with the experiment. However, if temperature 
goes lower further, it is predicted that Cp will exhibit a 
linear behavior at relatively lower temperatures, which is 
in need of experimental verifications. 
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APPENDIX A: FREQUENCY SUMMATION 

In this appendix, we shall first show that, when ^ > 0, 
the ideal Bose gas is still well defined within the extended 
ensemble theory, and then discuss the frequency sum- 
mation, which is very useful when we attempt to solve 
problems through Green's functions. 

Upon a transformation as in Eq. H3.2(l . the statistical 
average defined by Eq. (|2.3|) can always be transformed 
into a statistical average with respect to the Hamilto- 
nian H{b) of Eq. (|5.1(l . It is thus sufficient for us to 



show that the statistical average with respect to H{b) 
is properly defined when /i > 0. As is well known, ev- 
ery statistical average can be calculated through a cor- 
responding Green's function, and every Green's function 
can be derived from the generating functional of the sys- 
tem. Therefore, we need only to prove that the generat- 
ing functional is properly defined when ^ > 0. 

According to the path-integral formalism , the gen- 
erating functional W^[(/)+,0] for H{b) can be written as 
follows, 



J 



: J\f-^ J Vb+Vb expj i 



ihlOr, 



21.2 



2m 



5k,n - (pk.nO^ , 



(Al) 



where ujn — 2m:/ (3h [n S Z) denotes the Mutsubara frequency, and J\f the normalization factor. 



J\f=J I?6+2?5expji^6+„ 



(A2) 



If /I < 0, one gets immediately the familiar result. 



Wl 



cxp 



k,n 



\ 2m 



(A3) 



r 



Since fi,^k^/(2m) — /i > 0, the normalization factor J\f 
converges. 

When fi > 0, the normalization factor Af will diverge. 
However, that does no harm because the normalization 
factor can be cancelled by the numerator of Eq. lAlll . 
as is frequently encountered in quantum field theory [49j . 
After the cancellation, the rest is still Eq. (|A3() except 
/i > 0. As mentioned in Sec. [S] an unbounded function 
as integrand is permissible within Lebesgue integration. 
Therefore, the sum in the exponent of Eq. ljA3p is proper 
as a Lebesgue integral. The resulting ly 0] is also 
right because a function is permitted to be unbounded in 
Lebesgue integration. 

In a word, the generating functional for H{b) is always 
properly defined no matter how large the chemical po- 
tential /i is, which ends our proof. 

Now, there is no mathematical limit on the chemical 
potential of the ideal Bose gas, as is the case for the ideal 
Fermi gas. That is rational and significant, any limit on 
chemical potential should be set by the physical theory it- 
self rather than by the mathematical theory involved. In 
the extended ensemble theory, the chemical potential of 
a system is determined only by the physical requirement: 
the conservation of particles, whether the system obeys 
Bose-Einstein statistics or Fermi-Dirac statistics; there 



is no limit set by the mathematical tool, i.e., Lebesgue 
integration. 

As a consequence of Eq. (IA3|) , we obtain 



G (k, iuj„ 



1 



-2 X2 



t-1 //fkf. 
" V 2m 



(A4) 



where t/(k, iujn) represents the temperature Green's func- 
tion which is defined as [s^ 

g(k,r-r') = -TT{TMT)bl{T')}piHib))). (A5) 

As a function of k, Q(k,iLUn) can be unbounded if /i > 
0, but that will be all right because fj(k, ia;„) can only 
appear in the integrand of a Lebesgue integral over k. 

Those discussions indicate that, in the manipulations 
of +00 and — oo, one can benefit greatly from Lebesgue 
integration. That is because Lebesgue integration is de- 
signed, ab initio, on the set of extended real numbers: 
RtJ = Ru {-hoo} U {-oo} where M = (-oo, +oo) Ev- 
idently, those benefits can not be provided by Riemann 
integration. More benefits of Lebesgue integration will 
be seen from the following discussions. 
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When one evaluates the statistical average of an ob- 
servable, he will encounter the summation over Mutsub- 
ara frequencies. As an illustration, let us consider 



dk6[6k= / dkTY{blb^p{H{b))) 

JR3 



dk 



1 



Y.e'^-'^giKiLUn), (A6) 



where 77 = 0+. We shall discuss the thermodynamic limit 
case. Here, it is important to note that the frequency 
summation always goes ahead of the Lebesgue integral 
over k because the thermodynamic limit must be taken 
at last. 

To perform the frequency summation, it is helpful to 
make use of the complex function. 



1 



(A7) 



which has simple poles at the Mutsubara frequencies, i.e., 
z — iujn = i2mT / 0h, and transform the summation into 



a contour integral [s^ . 



1 



e^-"''g(k,*c.„) = J —e"^^f{z)g{k,z), 

(A8) 

where the contour C is depicted in Fig. 7. It should be 
emphasized that the above procedure requires that the 
function t/(k, z) must be analytic on the whole imaginary 
axis lest there be any poles other than iw„ on this axis. 
Sometimes, this requirement can not be met by every k, 
e-g-, 



V 2m 



M > 0. 



(A9) 



where E = IR^\D. Here the set D is removed form 
so that the function ^(k, z) can meet the requirement 
everywhere on E. 

Generally, ^(k, z) are analytic on both the upper and 
lower half complex-z planes [3^ Is^ . one can therefore 
deform the contour C into the contours F and C", which 
results in 

f i^e^..;(,)t;(k,,)^ f ii/(;,)c;(k,z). (AH) 
Jc 27rz Jc' ^ni 



z plane 



rZ = ICOn 





FIG. 7: The counters for the frequency summation of a Bose 
system. 



Each such Q(k,z) that satisfies |k| — ^/2mfi/h is singu- 
lar on z = 0, and thus can not meet the requirement. If 
such singular k's constitute just a null set D, i.e., C/(k, 2) 
satisfies the requirement almost everywhere on M.^ , which 
is the usual case, then those k's can be left out of con- 
sideration for they give merely a null contribution to the 
Lebesgue integral over k. Therefore, we can rewrite Eq. 
[a as 



= J^dkJ ^^e^^^f{z)g{k,z), 



(AlO) 



The contribution from the contour F vanishes owing to 
the convergence factor e*''^ [H^ . The right-hand side can 
be simplified as 



dz /■+°° 

f{z)g{k,z)=V dLufiLo)A{k,u;), (A12) 



where A(k., lo) is the spectral intensity. 



^(k, Lu)^~— [g{k, Lu + io+) - g{k, uj - io+)] 



and V represents the principal value. 



(A13) 



V 



+00 



du J{ijj)g{k,uj) 



lim 



diu J{uj)A{k,uj) 



+ 00 



d^/(c^)^(k,c^; 



(A14) 
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dk lim 



+0+ 



lim / dk 

a^0+ Je 



duj f{Lj)Aik,uj) 

y. 

duj f{uj)A{k,uj) 

dujfioj) [ dk^(k,t^) 
Je 

/ + 00 r 
dufiuj) / dk^(k,w). 



+ C30 



lim 

Q^0 + 



+00 



duj f{Lj)Aik,uj) 
duj f{Lj)Aik,uj) 

dLuficj) [ dkA{k,oj) 
Je 



(A15) 



Substituting Eqs. (jAllHA14() into Eq. HA10|) . one has 



Now, supplementing the set E with the null set one 
arrives at 

/+00 /• 
dLufiuj) / dkA{k,uj). (A16) 
-00 JR3 

As usual, by introducing the density of states A/'(w), 

AAH = -i-3 / dk^(k,c.), (A17) 
(27r) JR3 

Eq. (UTeli can be reduced into 




(A18) 



This is a very useful identity, as its applications, one can 
easily obtain the final results of Eqs. H5.35|l and H5.39|l . 
It should be stressed that Eq. (jAlSp holds whether the 
chemical potential /i is less than, equal to, or lager than 
zero. If At < 0, it reduces to the common result given 
in Ref. |54|. It is a generalization of the common result 
when fj, > 0. This generalization removes the mathemat- 
ical limit on the chemical potential. One can not reach 
the result of Eq. ljA18p within Riemann integration. 

Obviously, the above analyses for the ideal Bose gas are 
also suitable for other Bose systems, e.g., the interacting 
Bose gas, which will be studied in Sec. 17.11 

Remark: If the system obeys Fermi-Dirac statistics, 
the principal value is unnecessary. Instead of Eq. ljA18|l . 
one has 




(A19) 



I 

That is simply because the Fermi distribution function is 
bounded and continuous at w = for any finite temper- 
ature /3 > 0. A more direct reason for the unnecessity 
of the principal value consists in the fact that, unlike 
the Bose system, z = is not the pole of the complex 
function, 

because the Mutsubara frequencies for a Fermi system 
are all nonzero, i.e., cOn = {2n + 1) tt / fih 7^ (n G Z). 
As a consequence, the integration contour can goes as 
in Fig. 8, which gives straightforwardly the final result 
of Eq. HA19() . Evidently, that makes the Fermi system 
easier to handle mathematically than the Bose system. 
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FIG. 8: The counters for the frequency summation of a Fermi 
system. 
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APPENDIX B: LOW-TEMPERATURE 
THERMAL ACTIVATIONS OF A WEAKLY 
INTERACTING BOSE SYSTEM 

In Sec. 17.11 our discussions were confined within the 
Hartree-Fock approximation of Eqs. (|7.12l - l7!T5|l and the 
short-range approximation of Eq. H7.22|l . Within those 
approximations, one sees from Eqs. H7.40|l and (|7.44|) 
that the thermal activations at low temperatures can be 
expanded into even power series of T. We shall show that, 
beyond those approximations, this feature still holds as 
long as the interaction is weak and perturbative. 

Under the condition assumed, observables can be cal- 
culated via GF. The calculations can be transformed fi- 
nally into the integrals of the following form, 

J —oo ^ 



where g (lo^JT) is a real- valued function of w and /I, which 
comes from a summation of GF over k. As a function 
of u), g{uj,Jl) should be bounded at bottom or decrease 
fast with Lu so that the integral can converge. In both 
cases, we can cut off the integral from below at a certain 
frequency —lOc i^c > 0), that is. 



According to the definition of principal value, the right- 
hand side can be written as 



Observe that 



we have 



= -1, (B4) 



1 e-^*" - 1 



duj ^ ^ ' ^) _ _ f dujg(LU,u,)— [ duj ^ ^ ' ^\ + / do? ' . (B5) 

Integrated by substitution, it turns into 

9i^,'P) r% f ~, , rr. /•+°°, g{u;, + kBTx,]l) 
du—r, = — / doj q [LU — oJc, U] + KB-L / dx —3 

Jo ■ 

The second integral on the right-hand side vanishes exponentially as T — > 0, it can thus be neglected. The upper limit 
of the third integral can be set equal to +00 at low temperatures for the integrand decreases as when x becomes 
very large, 

, 9 r" , I ~- A 9 {ksTx, V)- g {-ksTx, Jl) 

"■^~bZ 7 = ^/ dujg{u; - uJcfJ-) + ksT dx . (B7) 

-00 — 1 Jo Jo e — 1 

The second integral on the right-hand side can expanded as 

> ^j{ksTx,l.)^^gi-ksTx,J^) ^ 2 g c(2n)g(^"-i) (0, ^) (fcBT)^""^ . (B8) 



Substituting it into Eq. IB 711 results in 



= - / g - c.,, m) + 2 ^ C(2n)5(^"-^^ (0, [ksTf^ , (B9) 
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the sum is an even power series of T . To the fourth order, it becomes 

^ d^fi^ZT = ^ [' dug{u;- Lo,, V-) + 2C(2)5' (0, Ji) {ksTf + 2C(4)g(3) (q, ^) (fc^T)^ , (BIO) 



it is sufficient for the discussion of the low-temperature 
properties of the system. This resuh impUes that the 
renormaUzed chemical potential Jl and the internal en- 
ergy E are even functions of T at low temperatures, 
which leads to the conclusion of Eq. H7.45(l : the spe- 
cific heat Cy of a weakly interacting Bose gas will vanish 
linearly as T ^ 0. 

The above result reminds us of the specific heat of 
the electron gas in a normal metal, which also vanishes 



linearly as T ^ 0. This linear behavior of the electron gas 
can be explained analogously. There, what is concerned 
is the integral. 



1' 



(Bll) 



where /i is the chemical potential. Along the same way 
as for the Bose gas, it can be expressed as 



/r "^""S^ " r " + 2 E (l - ^) C(2n)g(^"-i) (0, (fc^T)^" . (B12) 



To the second order, it is 



dLo^^^ - / dcogiu;- oj,, m) + ((2)^' (0, A^) (fc^T)^ , (B13) 



which is sufficient for the low-temperature properties of 
the electron gas in a normal metal because its Fermi tem- 
perature is very high, in comparison to the critical tem- 
perature of the BEC happening in a weakly interacting 
Bose gas. Eq. (|B13p also implies that the chemical po- 
tential and the internal energy E are even functions 
of T at low temperatures, and that the specific heat at 
constant volume Cy will vanish linearly as T — > 0. For 
example, one can easily verify, with the help of Eq. ljB13|l , 
that the /z, E and Cy for the ideal Fermi gas have the 
low-temperature forms, 



which are familiar results. Those discussions show that 
there exists similarity even between so different systems! 



/i(r) = /i(o) 
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